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ABSTRACT 


This research was designed to study the relationship between 
different amounts of vegetative ground cover and the energy reflected by 
corn canopies. Low altitude photography and an airborne multispectral 
scanner were used to measure this reflected energy. 

Field plots were laid out, representing four growth stages of corn. 
Two plot locations were chosen — on a very dark and a very light surface 
soil. Color and color infrared photographs were taken from a vertical 
distance of 10 m. Estimates of ground cover were made from these photo- 
graphs and were related to field measurements of leaf area index. Ground 
cover could be predicted from leaf area index measurements by a second 
order equation. 

Color infrared photography proved helpful in determining ground 
cover on dark soil backgrounds. Color photography w as useful for deter- 
mining ground cover on light soil backgrounds , as long as the ground 
cover did not exceed about T?%. 
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Microdensitometry and digitization of the three separated dye 
layers of color infrared film showed that the near infrared dye layer is 
most valuable in ground cover determinations. Computer analysis of the 
digitized photography provided an accurate method of determining percent 
ground cover. 

Multispectral scanner data were collected in two flights over the 
light soil background plots at an altitude of 305 m. Energy in eleven 
reflective wavelength bands from 0.U6 to 2.6 ym was recorded by the 
scanner. A set of eight ground reflectance panels was in close proxi- 
mity to the ground cover plots and was used to try to normalize the 
scanner data over time. 

Ratio techniques were used to relate unealibrated scanner response 
to leaf area index. The ratios of scanner data values for the 0.72 to 
0.92 ym band over the 0.6l to 0.70 urn band and the 1.0 to 1.1* ym band 
over the 0.6l to 0.70 ym band were calculated for each plot. The ratios 
related very well to leaf area index for a given scanner flight date, but 
could not be generalized between flights because of uncertainty in 
scanner response over time. 

Ground reflectance panels were used to relate laboratory reflec- 
tance measurements to scanner response. Separate prediction equations 
were obtained for both flight dates for all eleven reflective wavelength 
bands of the multispectral scanner. In this way, scanner response was 
normalized to ground panel reflectance. Ratios of normalized scanner 
data could be related to leaf area index over time. 
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The normalized scanner response was used to plot relative reflec- 
tance versus wavelength for the ground cover plots. Spectral response 
curves resulted which were similar to those for hare soil and green 
vegetation as determined by laboratoiy measurements. The spectral res- 
ponse of different ground cover plots represented a "mixing" of the 
spectral response curves for the bare soil and green vegetation compo- 
nents of the scene. 

The spectral response curves from the normalized scanner data 
indicated that reflectance in the 0.72 to 1.3 vm wavelength range in- 
creased as leaf area index increased. A decrease in reflectance was 
observed in the O.65 um chlorophyll absorption band as leaf area index 
increased. This confirmed the validity of using the ratio of the res- 
ponse from a near infrared vavelength band to that of the red wavelength 
band in relating multispectral scanner data to leaf area index in corn. 
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INTRODUCTION 

Properties of reflectance and emittance of energy have been studied 
for many natural earth surface features in connection vith research in 
the field of remote sensing. Usually the spectral properties of indivi- 
dual samples of soils or plant materials are determined in the labor- 
atory* and their properties axe generalized to the field situation. Very 
little work has been done in measuring the spectral properties of the 
green vegetation - bare soil complex in the field. 

An understanding of the component spectral properties of bare soil 
and green vegetation is a first step in determining the spectral pro- 
perties of a crop canopy. Incident solar spectral properties also must 
be taken into account in the field situation. 

In this study* photography and a multispectral scanner were used 
to detemine some of the spatial* spectral, and temporal aspects of in- 
creasing ground cover in com crop canopies. The com canopy was chosen 
because of the prevalence of the com crop in the Midwestern United 
States , and because of the interesting geometrical complexity of the 
canopy. Percent ground cover has always been most difficult to measure 
in com canopies and has been a measurement of questionable value to the 
agronomist. Leaf area index could also be considered a measurement of 
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vegetative cover, and is of great value to the agronomist in determining 
the extent of utilization of energy from the sun in photosynthesis. 

The potential for measuring vegetative cover from the spectral 
properties of a canopy holds promise in determining the above ground 
biomass of a crop by remote sensing techniques . In a general way the 
above ground biomass could then be related to crop yield. 

Future efforts in remote sensing from orbital altitudes such as are 
proposed for the Earth Resources Technology Satellite (ERTS) and SKY LAB 
will be concerned with general views of agri cultured crops. With the 
extremely high altitude and coarse resolution from space platfoms , it is 
likely that any discrimination of healthy green agricultural crops will 
be on the basis of differences in vegetative cover, and not on in- 
dividual plant spectral properties. Vegetative cover may well reflect 
disease infestations , drought stress , and other physiological stresses 
within a given crop type. 



CHAPTER I 


REVIEW OP LITERATURE 


Crop Canopy Considerations In Remote Sensing 

Remote sensing has been developing in recent years into a useful 
tool for the agronomist in characterizing the spectral, spatial, and 
temporal aspects of soils and plant communities (Hoffer, 1967; LARS, 
1967; LARS, 1968; LARS, 1970). The extent of vegetative cover has been 
an important consideration to the crop and soil scientist in viewing 
agricultural land from a remote position. Multispectral scanner data 
from fields of bare, uniformly cultivated soil have been useful in de- 
lineating soil boundaries of interest to the soil mapper. Scanner data 
from full crop covers have been useful in general crop inventory and 
detection of physiological stress in plants. A partial ground cover 
tends to complicate the interpretation of the scene for both the crop 
and soil scientist. Knowledge of the extent of vegetative cover will be 
necessary as agronomists attempt to refine their interpretation of re- 
motely sensed data. 

Remote sensing has potential utility in the determination of the 
amount of vegetative cover in the above ground portion of the crop bio- 
mass. Miller and Pearson (1971) attempted to measure spectrally the 
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standing crop biomass of a shortgrass prairie vegetation. They reasoned 
that a measure of percent vegetative cover or leaf area index would 
relate to the carrying capacity for grazing in an area. It is reasonable 
to assume that the same principle could be used to relate vegetative 
cover in corn canopies to the above ground biomass of corn and, perhaps, 
in a relative way to corn yield. 

Spectral response data obtained with an airborne multispectral 
scanner over a particular crop represent an integration of the radiation 

from both plants and soils in the instantaneous field of view of the 

» 

scanner. Many variables come into play which affect the Intensity and 
wavelength of radiation from both plants and soil. Leaf angle, leaf 
moisture content, leaf area index, plant geometry, variety or species, 
and planting technique will affect the spectral response from plants in 
a field situation. Moisture content, texture, surface crusting, organic 
matter content, surface residue, and color will affect the spectral 
response from surface soil. In order to understand the soil-crop complex 
it is first necessary to look at the individual factors which have an 
effect on intensity and wavelength of the energy radiating from plants 
and soils. 


Plant Reflectance Properties 

The spectral response from an individual plant leaf in the O.U to 
2.6 urn wavelength reflective range can be broken up into three portions 
(Hoffer and Johannsen, 1969): 1, The visible wavelengths where pi©nent 
absorption predominates, 2. The 0.72 to 1.3 vm wavelength region where 
very little absorption takes place and almost all the incident light is 
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transmitted or reflected, and 3. The 1.3 to 3.0 ym wavelength region 
where prominent water absorption takes place at 1.45 ym and 1.95 ym. 
Chlorophyll, carotene, anthoeyanin, and xanthophyll absorb at 0.445 ym 
in the blue, while chlorophyll alone absorbs at 0.645 ym in the red 
portion of the spectrum (Gates, et al. , 1965). Reflectance in the vi- 
sible peaks at 0.540 ym in the green region for green leaves and accounts 
for the green color of plants perceived by the human eye. According to 
Khipling (1970) about 40 to 60 percent of the infrared radiation in the 
wavelength range from 0.7 to 1.3 ym is scattered upward through the sur- 
face of incidence and is designated as reflected radiation , whereas the 
remainder is scattered downward and transmitted. The strong absorption 
by a leaf in the infrared beyond 1.3 ym is due to water » as supported by 
theoretical considerations of the equivalent water thickness of leaves 
(Knipling, 1970). Absorption of radiation in plants is very high in the 
visible portion of the spectrum where energy is needed for photosynthesis 
and drops quite low from 0.7 to 1.1 ym where most of the energy of the 
solar spectrum is concent rated. Absorption rises again in the 1.3 to 
3.0 ym region. Thus, a plant's spectral properties allow it to absorb 
energy where it is needed for growth, and reflect and transmit energy 
where it is detrimental. 

Leaf moisture content hem a consistent effect on reflectance in 
all reflective wavelength regions. Johannsen (1969). working with corn 
and soybeans , found a high negative correlation between DK-2 spectro- 
photometric leaf reflectance and leaf moisture in the chlorophyll absorb- 
tion band at 0.64 ym, and concluded that leaf moisture may be changing 
leaf pigments in short periods of time. Johannsen also found high 



negative correlation betveen leaf reflectance and leaf moisture in the 
water absorption bands at 1.43 Pm and 1.94 pm and in the green peak at 
0.53 um. He concluded that reflectance measurements from corn and soy- 
bean leaves showed a steady increase in all reflective wavelengths with 
decreasing leaf moisture. Johannsen reported diurnal effects on leaf 
moisture content. Leaf moisture seemed to increase with increasing soil 
moisture. Scanner data from com canopies gave decreasing reflectance 
with increasing plant and soil moisture. 

Leaf structural parts exhibit varying effects on reflectance pro- 
perties of plants. Gausman, Allen, and Cardenas (1969), working with 
salinity stressed cotton, reported that changes in reflectance in the 
0.75 to 1.3 pm range were primarily due to changes in internal structure 
and not leaf moisture content. They associated increased reflectance 
with thicker leaves having stronger palisade cell development , loosely 
arranged spongy mesophyll, and hence, more intercellular space or inter- 
faces. Sinclair (1968) found that the greaitest reflectance in the 
shorter infrared wavelengths came from the internal leaf structure which 
had the greatest amount of cell wall material perpendicular to the inci- 
dent radiation. Thus, he found that the palisade tissue was more reflec 
tive in dorsiventral leaves (typical of dicotyledonous plants like the 
soybean) since it provided greater barriers of cell walls to impinging 
radiation than did spongy mesophyll tissue. Sinclair also reported that 
the cuticle probably contributes a relatively small portion of the total 
leaf reflectance observed in the wavelength bands highly reflected by 
leaves. In the highly absorbed bands, though, he found that cuticular 
reflectance may contribute nearly all of the leaf reflectance observed. 
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Plant leaves may be diffuse or specular reflectors depending on the 
wavelength of interest. Sinclair (1968) reasoned that in order to obtain 
high levels of reflectance in the shorter infrared wavelengths, the cell 
walls must be diffuse reflectors of incident radiation. Breece and 
Holmes (l97l)» in a study of the bidirectional scattering characteristics 
of corn and soybean leaves, found that for wavelengths beyond 0.75 pm, 
reflectance was nearly Lambertian in character. They reported that the 
more strongly absorbed visible wavelengths showed a tendency to reflect 
specularly the radiation incident upon the leaf. 

When the spectral response from stacked leaves or leaves in a 
canopy is considered, the idealized single leaf spectral reflectance 
becomes much more complicated. Jfyers et al. (1966), in an experiment 
with stacked leaves , reported that no difference in reflectance occurred 
in the visible wavelengths for any combination of stacked leaves , but a 
17 % increase in the reflectance occurred in the near infrared for two 
leaf layers. As the number of leaf layers increased up to six layers, 
the reflectance increased in the near infrared except in the 1.45 pm and 
1.95 Pm water absorption bands. The probable explanation for this en- 
hanced reflectance with multiple leaf layers in the near infrared is a 
retransmission of infrared radiation back through the canopy from lower 
layer reflection. Jfyers (1970) reasoned that the leaf layers could be 
related to an actual measure of leaf density or leaf area index (LAI). 
Leaf area index is the ratio of leaf area to soil area in a plant canopy. 
Jfyers reported that Kodak Ektachrome Infrared Aero film, which is sensi- 
tive in the wavelength interval from 0.4 to 0.9 pm, could be used to 
record variations in reflectance associated with leaf area index. 
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For the real situation of leaves in a crop canopy, reflectance is 
considerably less than for a single leaf because of general attenuation 
of radiation by variations in illumination angle, leaf orientation, 
shadows, and non foliage background (soil in most cases). Vinogradov 
(1969) concludes that the reflectance for vegetation decreases with an 
increase in the roughness of the surface geometry of the plant canopy. 

He states that the reflectivity for the soil-plant system decreases with 
the increase in leaf area index, but that beyond an LAI of 3 or 4, a 
further decrease in the reflectance practically does not occur. 
Vinogradov, however, reports that reflectance in the 0.72 to 0.98 um 
wavelength band is augmented with an increase in ground cover, height, 
and leaf area index. Knipling (1969) estimates that the visible and 
near infrared reflectance from a complete canopy may be 3# to 5# and 35# 
respectively versus 10# and 50# for a single leaf. This would be a 60# 
reduction in reflectance in the visible and a 30# reduction in reflec- 
tance in the near infrared. The relatively smaller reduction in the 
near infrared is accounted for by enhanced reflectance in the multiple 
leaf layered-canopy. 

Reflective Properties of Background — 

Soil , Shadows , and Hon-Green Vegetation 
Baumgardner et 8 lL. (1970) found that organic matter plays an 
important role in bestowing spectral properties to soils when the organic 
matter content is above 2#. As the organic matter content drops below 2# 
it may be less effective in maisking out the effects of other soil consti- 
tuents such an iron oxides. 
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Boiyoucos (1913) showed that different types of soil tended to 
radiate heat energy differently when dry and about the same or more when 
well moistened and in their natural condition. Hoffer and Johannsen 
(1969) showed that moist soils had an overall lower reflectance than 
their dry counterpart in the O.U to 2.6 um wavelength region. They also 
found, however, that crusting of surface soil earn make a soil appear dry 
when it is actually wet. Cipra et al. (1971) found that crusted surfaces 
gave higher reflectance values in the O.U 3 to 0.73 um wavelength region 
than did soils with the crust broken. They attributed the lower reflec- 
tance of the disturbed soil to the rough surface which presumably caused 
scattering of light as well as a shadowing effect. Because of this 
crusting effect the emissive wavelengths in combination with the reflec- 
tive wavelengths could give the best clues to soil moisture content. 
Bowers and Hanks ( 1965 ) estimated that an additional l4# to 16% of the 
direct solar energy was absorbed by a silt loam due to an increase in 
moisture content. 

Differences in soil texture may account for differences in spectral 
response. Bowers and Hanks ( 1965 ) found differences in reflectance as 
large as 2k% evident for soil particles ranging from 22 um to 2680 um in 
size (coarse clay to sand size fractions). They found a rapid exponen- 
tial increase in reflectance with decreasing particle size. The field 
situation may not be very similar since surface soil aggregates of vary- 
ing sizes will predominate. 

Soil color may be related to organic matter content, the presence 
of oxidized or reduced iron compounds, and the internal drainage char- 
acteristics of the soil. Lighter colored soils may have lower organic 
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matter contents, higher percentages of oxidized iron compounds, and good 
natural drainage. Darker soils may have high organic matter contents, 
reduced iron compounds, and poor natural drainage. For the 0.3 to 1.0 um 
wavelength region, Condit (1970) studied the spectral properties of many 
surface soils having a wide variety of colors and textures. 

Factors whose effects on the spectral reflectance of surface soils 
axe most important, yet more difficult to quantitize, sure soil surface 
roughness , shadowing, and non-green crop residue. Tillage effects and the 
resulting clod size of cultivated soil will impart varying spectral 
characteristics to bare soil. As for shadow effects, Miller and Pearson 
(1971) felt that shadows might be considered a spectrally separate 
material upon first investigation , but that shadows would more accurately 
be described in terms of modified spectral properties of the materials 
upon which the shadows fall. Shadows in a canopy are certainly elusive 
and hard to measure at ground level. Non-green crop residue will have 
to be considered as part of the background scene when it is present in 
large enough amounts to affect the spectral properties of the canopy. 
Although this component may play a large part in grassland canopy con- 
siderations , non-green vegetation should not be too complicating a factor 
in cultivated crops. 

Incident. Solar Spectral Properties 

The character of reflected and emitted light from a plant canopy 
will depend on the nature of the surface and the intensity and spectral 
quality of the radiation incident upon the surface. The plant canopy may 
be irradiated by direct sunlight, scattered skylight, and thermal radia- 
tion from the surroundings. Idso, Baker, and Gates (1966) reported that 
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scattered or diffuse radiation (skylight) makes up as much as 10 % to 15 % 
of the direct sunlight on a clear day. Cloud light, or the diffuse 
radiation which is the result of a complete overcast, is another type of 
solar radiation which may be present. Idso (1966) also lists two other 
types of solar energy: 1. that transmitted through the leaves of a 

canopy, and 2. that radiation reflected from other plants (which is 
usually rich in the green and near infrared wavelengths). Figure 1 from 
Gates (1962, 1963, 1965) shows the distribution of solar energy at the 
earth's surface with the accompanying absorption bands for 0 2 » 0^, and 
HgO. Figure la. is in the familiar increments of wavelength while 
Figure lb. is incremented in wave numbers. The various absorption bands 
are evident, and the great variation between the intensity and quality of 
direct sunlight, cloud light, and skylight is made evident. The choice 
of spectral bands for analysis of a vegetated surface would have to be 
made with reference to the incident ^ solar energy as well as the soil and 
plant spectral characteristics. 

Methods of Estimating Percent Vegetative Cover 
Direct Measurements 

According to Miller and Pearson (1971) there is no rapid method 
which has been found practical for measuring percent ground cover from 
the ground location of the observer. Estimates of percent ground cover 
can be made by individual research workers and checked for consistency 
with estimates of other researchers. The subjectivity involved in ground 
level estimation of percent ground cover is even further increased in a 
complex canopy such as the corn canopy where view angle effects, height 




WAVELENGTH, MICROMETERS 



b. wave number increments 


Figure 1. Solar energy spectrum; incremented in wavelength units and 
wave numbers (Gates, 1962, 1963, 1965). 
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effects, and shadowing complicate the scene. Nevertheless, ground 
measurement methods for percent ground cover have been attempted, and 
Cuellar (1971) suggests methods for two different canopy structures: 


1. Solid canopy with bare soil exposed only in the 

furrows: The "bare soil width" and row width are 

measured, the bare soil width being the amount of 
bare soil showing between the canopies of adjacent 
rows . Percent crop cover is calculated as : 

, (RW-BW) x 

— — ) X 100 = percent cover 

where 

FW * row width in cm. 

BW a bare soil width in cm. 

2. Open canopy with bare soil exposed through the 

canopy and furrow: The open canopy is considered 

as being "solid," and the formula above is used. 

A subjective estimate is made of the open spaces 
in the canopy, and this amount is subtracted from 
the assumed cover to obtain actual cover. 


In the second method, used for corn crop canopies, the subjectivity 
involved can be readily seen. Cuellar found it especially difficult to 
estimate cover as ground cover approached 100£, especially in tall crop 
canopies such as corn where ground observers were unable to look down at 
the crop. 

A more quantitative measure for characterizing the extent of vege- 
tative cover would be that of leaf area index (LAI). Although leaf area 
index has not yet been found to relate directly to percent ground cover, 
its usage is consistent with multi-layer models of plant canopies , as 
discussed earlier. Leaf area index, of itself, should certainly be a 
valuable measure of above ground biomass . In the case of the corn 


canopy, the LAI is easily calculated by measuring the length and width of 
each leaf of a corn plant , multiplying these together times a factor of 
0.73 to account for the non-rectangular form of the com leaf (McKee, 

1 964), and multiplying this times the plant population per unit area of 
soil. 

For shortgrass prairie vegetation the point -quadrat method has 
proved valuable for measuring percent ground cover (Miller and Pearson, 
1971). This method would be impractical for canopies of considerable 
height. Equally unsuitable would be destructive sampling methods of the 
canopy vegetation when additional measurements are desired over time. 

Photographic Methods 

Various photographic approaches have been used to arrive at an 
estimate of the amount of surface covered by a vegetative canopy. 

Perhaps the most unusual approach was that of McCree (1968) who chose to 
view the plant canopy from the ground level looking vertically upwards. 

He used a 180° mirror to photograph the canopy of willow trees and was 
able to relate tonal response on color infrared film with sunlit leaves , 
clear sky, shade light, and sunlight. Similar photographs have been 
taken of com canopies (LARS, 1968) from a ground position with a wide 
angle lens. The resultant photographs give a qualitative view of the 
canopy cover and the amount of transmitted light. 

In viewing the plant canopy from a vertical position, Miller (1969) 
found that dark colored soils and dark, wet soils had a suitably low re- 
flectance in the 0.7 to 0.9 pm wavelength band so that response in the 
instantaneous field of view of a crop canopy would be approximately equal 
to percent cover of vegetation. Hoffer (1967) also found the 
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photographic in frared wavelengths useful in estimating ground cover on 
dark soils and noted that increasing canopy resulted in higher response 
on infrared film. For light soils , with high reflectance relative to 
plants, the visible photographic wavelengths were best for determining 
percent vegetative cover. On light soils the 0.7 to 0.9 vm band showed 
a blending of soil and highly reflective vegetation resulting in a 
relatively unifora response despite variations in ground cover. In the 
visible wavelengths, however, the light colored soil did not blend with 
the vegetation and the density of vegetative cover was proportional to 
the decrease in response registered on film. 

A sampling method was used by Null (1969) to measure canopy density 
from photography. He used a point grid system to estimate the percent of 
film area represented by vegetation. 

Spectrum Matching 

For estimating percent green vegetation on light colored soils , 
Miller (1969) suggests a technique of "spectrum matching" in which pure 
samples of vegetation and soil are "mixed" to fora a set of test spectra 
representing various levels of percent vegetative cover. The spectral 
response values from the instantaneous field of view of the remote 
sensing device could be matched with these test spectra to estimate 
percent ground cover. Figure 2 (Miller, 1969) illustrates this proce- 
dure for low, medium, and high percent ground cover situations. 
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Ratio Techniques 

Kristof and Baumgardner (1970) found that a decrease in radiation 
in the visible and an increase in the reflective infrared region occurred 
in comparing multispectral scanner data over crop cover from April to 
late June. Disregarding complications of row structure and direction 
which may produce a "stove pipe effect" (LARS, 1968) in corn canopies, 
the increase in ground cover could be detected by establishing a ratio 
of the relative reflectance in the chlorophyll absorption band to the 
relative reflectance in the near infrared. Kristof and Baiungardner 
(1970) used the ratio of relative reflectance in the O.58 to 0.62 urn band 
to the relative reflectance in the 0.80 to 1.00 urn band to observe in- 
creases in ground cover. With the increased reflectance in the near 
infrared and the increased chlorophyll absorption with increasing crop 
cover, this ratio would decrease with increasing crop canopy. For a 
given crop and soil type an estimate of the percent ground cover could 
possibly be made from the value of the ratio of the relative reflectance 
from two channels of a multispectral scanner. 

Theoretical Formulations 

Various investigators have attempted theoretical models for esti- 
mating the radiation response from the soil background and plant com- 
ponents of a plant canopy. Allen and Richardson (1968) used the two 
parameter Kubelka-Munk theory to arrive at theoretical values for re- 
flectance and transmittance of a plant canopy from the leaf area, an 
absorption coefficient, a scattering coefficient, and the background 
reflectivity. The two coefficients are based on crop geometry and the 
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optical properties of the individual leaves. The validity of this theory 
in an actual situation is doubtful, since it assumes a uniform distri- 
bution of leaf inclinations and assumes infinite lateral extension to the 
canopy. 

Idso and deWit ( 1970) improved on the Kubelka-Munk theory by 
including a factor called the leaf distribution function, which char- 
acterizes the frequency distribution of plant leaf inclination from the 
horizontal. The so-called deWit-Idso theory introduces leaf area index 
and canopy density to consider the extinction of the solar beam as it 
penetrates the canopy foliage. The deWit-Idso theory gives information 
on the fraction of leaves in each layer receiving diffuse and direct 
light at various angles of inclination to the sun's radiation. 

Gates (1965) arrived at a formula for the total radiation absorbed 
by a stand of vegetation and soil. Figure 3 from Gates (1965) shows that 
the energy absorbed by a plant canopy increases with exposed soil up to 
a point (about ho% ground cover) but then declines in a parabolic manner. 
Gates found that the reflected solar component has only a fairly strong 
dependence upon the amount of plant cover on dark soils. 

Holmes and MacDonald (1969) arrived at a formula for the integrated 
radiance vector for crop canopies as follows : 

total radiation = 0 +6 +Y +6 

r BL r SL r BS r SS 

where r = the radiation vector for: BL (bright leaves) 

SL (shaded leaves) 

BS (bright soil) 

SS (shaded soil) 

a »8»Y>6 are terms for the ratios of the projected areas of the four 

component scenes to the projected viewed scene area, on a plane 

normal to the line-of-sight. 



Q„bs = 068 +0.42f -0361 



Figure 3 * Energy absorbed by a plant canopy as a function of exposed soil (Gates, 1965) 
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The sum of the four component ratios should be one. This formula is 
valuable for its inclusion of the components of shaded leaves and soil 
in the radiation scheme. Variations in view angle, sun angle, and row 
direction are considered in this formulation of crop canopy radiance as 
detected by the instantaneous field of view of an airborne scanner. 

A very thorough review of plant and soil reflective properties was 
conducted by Kumar (1972). This review contains a more detailed des- 
cription of theoretical plant canopy models. Krumpe (1972) compiled a 
bibliography of over 850 references on remote sensing of terrestrial 
vegetation, covering most of the material published from 1955 to the 
present . 
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CHAPTER II 

MATERIALS AND METHODS 
Plot Location 

During the summer of 1971 an experiment vas carried out utilizing 
the facilities of the Laboratory for Applications of Remote Sensing 
(LARS) and the Purdue University Agronomy Farm to study the extent of 
vegetative cover of different growth stages in com canopies. The ob- 
jective was to observe and study the reflectance characteristics of com 
canopies with both light and dark surface soil backgrounds. To this end, 
two plot locations were chosen at the Purdue University Agronomy Farm 
representing a light colored forest soil and dark colored prairie soil. 
Agronomy Farm plot number 15 was chosen as the site with a light surface 
soil, Russell silt loam, a well-drained alfisol with a dry Munsell color 
of 10 YR 6/2 and an organic matter content of about 3 %> Agronomy Farm 
plot number 44 was chosen as the site with a dark surface soil, Chalmers 
silty clay loam, a very poorly drained mollisol with a dry Munsell color 
of 10 YR 4/1 and an organic matter content of about 6 %. 

The plots were located, on the northern edges of plots 15 and 44, 
both bulk com plots, so that shadowing problems would be minimized when 
data were taken from a Hi-Ranger bucket above the plots. The light 
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colored Russell soil site was very near the Agronomy Farm Weather Station, 
location of LARS ground reflectance panel deployment, and was thus in the 
path of several flights by the University of Michigan's multispectral 
scanner. The scanner plane did not fly over the dark colored Chalmers 
soil site. 


Plot Design 

In order to facilitate data collection and provide a wide range of 
ground covers on any given date of a scanner overflight or photographic 
mission, four different planting dates were used. Pioneer 3369A corn was 
planted in three replications at each site location bn each of the 
following four dates: April 30, May 22, June 9, and June 30. This pro- 

vided a spacing of 22, 18, and 21 days between planting dates, respec- 
tively. 

Each plot extended 55 meters along the edge of the bulk corn fields 
and consisted of twelve sub-plots 1*. 57 by 4.57 meters in size. Each sub- 
plot contained six rows of corn, 4.57 meters long in 7 6 cm rows. Repli- 
cations of planting dates were assigned randomly to the twelve sub-plots. 
Table 1 illustrates the layout of the two plots. 

No special fertility treatments were used, and adequate fertility 
levels were maintained for the Agronomy Farm's bulk corn production. A 
plant population of about 49,400 plants per hectare was desired, and 
actual populations of from 43,200 to 50,700 plants per hectare were ob- 
tained. Act vial plant counts were taken in each sub-plot for the calcu- 
lation of leaf area within the sub-plot. 


Table 1. Assignment of planting dates to field plots. 

Planting Dates 

A = April 30 
B - May 22 
C = June 9 
D» June 30 
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Data Collection 

At the time of scanner overflights and photographic missions, 
field measurements of plant height and leaf area were made. Two plants 
were chosen at random from each sub-plot , and the plant height and leaf 
area were calculated from the average of these measurements. Leaf area 
was calculated using the method of McKee (1964) in which the length and 
width of each leaf on the corn plant are measured and multiplied together 
times a factor of 0.73 to account for the non-re ct angular form of the 
com leaf. Leaf area index is then calculated from the average leaf area 
per plant and the plant population per unit area of soil. 

On the four dates, July 13, July 21, August 3, and August 17, the 
LARS Hi-Ranger truck was used to obtain vertical photographs over the 
sub-plots from an altitude of approximately 10 meters. Color and color 
infrared photography were taken on these dates with two Nikkormat 35 mm 
cameras. A typical frame of photography contained four full rows of each 
sub-plot, so that border effects could be eliminated. 

Numerous multispectral scanner overflights were available from the 
summer of 1971* but only two of these were selected for analysis. The 
University of Michigan's multispectral scanner collected spectral data in 
eleven channels in the reflective wavelength range of 0.46 to 2.60 urn 
from overflights on July 12 and July 21, 1971. Both flights were at an 
altitude of about 305 meters and were flown in a north-south flight path 
directly over the Agronony Farm Weather Station. The scanner data tapes 
were converted from analog to digital form and were reformatted for pro- 
cessing on the LARS IBM 360/67 digital computer. 
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CHAPTER III 

LOW ALTITUDE PHOTOGRAPHY OF CORN CANOPIES 
Problem Definition 

Numerous investigators have examined the usefulness of color and 
color infrared photography of plant canopies in determining percent 
ground cover (Hoffer, 1967; McCree, 1968; Miller, 1969). Color film is 
generally felt to be best for ground cover determination on light colored 
soil backgrounds. Color infrared film has been of value in ground cover 
determination on dark colored soil backgrounds. 

Characteristics of color and color infrared film prevent the re- 
searcher from obtaining quantifiable spectral information about plant 
canopies from these film types. Relative spectral measurements can be 
made, however, by photographic separation of the dye layers of the film 
and microdenaitometry of these color separations. Computer analysis of 
digitized photography provides a means for distinguishing between the 
components of green vegetation and bare soil by analysis of spectrally 
separable classes. Digitization of low altitude photography thus pro- 
vides a useful technique for analysis of the spectral and spatial re- 
lationships in plant canopies. 
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Procedures 

Photographic Data Collection and Ground Measurements 
Vertical photographs were obtained from the 2h individual sub-plots 
on four dates: July 13, July 21, August 3, and August 17, 1971. The 

LABS Hi-Ranger truck was used to take the photographs from an altitude of 
about 10 meters. Color (Kodak Ektachrome-X and Kodachrome II) and color 
infrared (Kodak Ektachrome Infrared Type 2^3) transparency films were 
used in two Nikkormat 35 mm cameras with 50 mm focal length lenses to 
obtain the photographic data. A Kodak Wratten No. 12 filter was used 
with the Kodak Ektachrome Infrared film, and a skylight filter was used 
with the Kodak Ektachrame-X and Kodachrome II films. Simultaneous 
photographs were taken of each plot with the two cameras in such a way 
that four com rows filled the field of view. The resulting scale on 
the 35 mm transparencies was 1:200. 

Within one day of the photographic data collection, ground measure- 
ments were made of plant height and leaf area. Two plants were chosen 
at random from the sub-plots for these measurements , and leaf area was 
determined by the method described in the previous chapter. No attempt 
was made to measure actual percent ground cover since the geometry of the 
canopy prohibited anything but rough estimation of ground cover from 
ground level. Leaf area index determinations were preferable because of 
the ease and non-subjectivity of measurement. Percent ground cover es- 
timates were limited to measurements taken from the low-altitude photo- 
graphy using the dot grid system described in the next section. 
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Estimation of Percent Ground Cover by Point Grid System 
Thirty photographs were selected for percent ground cover esti- 
mation by a point grid technique. This included 12 Kodachrome II photo- 
graphs of the Russell plots taken July 13, 12 Ektachrorae-X photographs 
of the Russell plots taken July 21, 3 Ektachrome-X photographs of the 
Russell plots taken August 3, and 3 Ektachrome Infrared photographs of 
the Chalmers plots taken on July 13. These particular photographs were 
chosen because of their high quality and because of the diversity of 
ground covers represented among them. 

The 2k by 3 6 mm transparencies were projected onto a rear-viewing 
screen and enlarged 15 times to a size of 1 6 by 2k cm. This gave a 
scale of about 1:30 for the enlarged image. A point grid of 1 point per 
5.6 mm was overlayed on the screen, and a rectangular sample of two corn 
rows from the center of the frame was outlined. All the points which 
fell on com plants were tallied and expressed as a percent of the total 
number of points within the bounded area, resulting in an estimate of 
percent ground cover. Table 2 lists the 30 plots , their planting dates , 
observation dates, leaf area indices, and point grid estimates of percent 
ground cover. Regression analysis was performed to 'determine if percent 
ground cover could be related to leaf area index. 

Analysis of Digitized Photography 
Characteristics of Color Infrared Film . The range of spectral 
sensitivity of color infrared film is extended into the near infrared 
wavelength region, where vegetation typically has a very strong spectral 
response. Kodak Ektachrome Infrared Film Type 2kk3 t is the latest 
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Table 2. Point grid estimates and leaf area index measurements for 
30 ground cover plots. 


Estimated 


Plot 

Planting Date 

Observation Date 

LAI 

Percent Ground Cover 

RP1 

June 30 

July 13 

0.01 

3.3 

RP2 

May 22 

July 13 

2.89 

87.2 

RP3 

April 30 

July 13 

3.48 

93.8 

RP4 

June 9 

July 13 

1.30 

6o.o 

RP5 

June 30 

July 13 

0.01 

3.3 

RP6 

May 22 

July 13 

2.34 

88.8 

RP7 

June 9 

July 13 

1.11 

52.8 

RP8 

June 30 

July 13 

0.04 

7.6 

RP9 

April 30 

July 13 

3.10 

87.2 

RP10 

June 9 

July 13 

1.17 

42.8 

RP11 

May 22 

July 13 

2.93 

83.2 

RP12 

April 30 

July 13 

3.81 

91.8 

RP1 

June 30 

July 21 

0.11 

3.5 

RP2 

May 22 

July 21 

3.00 

94.6 

RP3 

April 30 

July 21 

4.28 

89.5 

RP4 

June 9 

July 21 

2.31 

81.5 

RP5 

June 30 

July 21 

0.12 

4.8 

RP6 

May 22 

July 21 

3.57 

89.1 

RP7 

June 9 

July 21 

2.15 

78.2 

RP8 

June 30 

July 21 

0.17 

12.2 

RP9 

April 30 

July 21 

3.38 

93.9 

RP10 

June 9 

July 21 

2.10 

72.7 

RP11 

May 22 

July 21 

3.45 

88.5 

RP12 

April 30 

July 21 

4.06 

94.7 

RP1 

June 30 

August 3 

0.57 

29.0 

RP5 

June 30 

August 3 

0.63 

23.8 

RP8 

June 30 

August 3 

0.93 

34.8 

CP9 

April 30 

July 13 

4.44 

96.5 

CPU 

June 9 

July 13 

1.31 

70.7 

CP12 

May 22 

July 13 

2.93 

90.3 


RP = Russell Plot 

CP ** Chalmers 

Plot 
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version of Kodak's multi-emulsion color infrared film. Fritz (1971) 
illustrates the wavelengths of spectral sensitivity for the film's three 
emulsion layers as well as the cutoff point for the Kodak Wratten No. 12 
filter (Figure U). 

Knipling (1969) describes the spectral sensitivities of color 
infrared film and the image forming processes. The film has three dye 
layers: 1. a cyan forming layer sensitive to the near infrared wave- 

lengths* 2. a magenta forming layer sensitive to the red visible wave- 
lengths, and 3. a yellow forming layer sensitive to the green visible 
wavelengths. All three emulsions are sensitive to radiation below 0.5 um, 
but these wavelengths are eliminated using the deep yellow Wratten No. 12 
filter. Thus, color infrared film has two dye layers sensitive to 
visible light excluding the blue, and one dye layer sensitive to the near 
infrared. It is stressed that this film is not sensitive to the thermal, 
but rather is sensitive only up to 0.89 um in the shorter wavelength 
infrared. Ektachrome Infrared Film 'type 2M3 is sensitive to wavelengths 
from 0.1*7 to O.89 um as compared to Ektachrome-X or Kodachrome II which 
are sensitive to wavelengths from 0.1*0 to 0.71 um. The 0.72 um wave- 
length is generally considered to be the boundary between the visible 
and infrared wavelength regions. 

Color Separation . Twenty color and color infrared transparencies 
were chosen for color separation and microdensitometry. Sixteen color 
and four color infrared transparencies were among those used. Before the 
color and color infrared multi-emulsion films could be densitometered, 
they had to be separated by a standard commercial process into the three 
emulsion layers and reproduced as black and white positive transparencies. 



Sensitivity Transmittance (%) 
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0.4 0.5 0.6 0.7 0.8 0.9 

Wavelength (y Ltm ) 

Figure U. Spectral sensitivity curves for Kodak Ektachrome Infrared 
film type 24^3, and the spectral transmittance of Kodak 
Wratten filter No. 12 (Fritz, 1971). 
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Hoffer, Anuta, and Phillips (1971) describe this separation process as 
well as all the other steps in digitization of photography. Kodak 
Wratten No. l*7b blue, No, 6l green, and No. 23a red filters were used 
for the separations. Eastman Kodak Pamphlet No. E-59 (1968) describes 
the use of Wratten filters for emulsion layer separation of color films. 
Table 3 lists the bands of spectral sensitivity of the black and white 
transpanencies upon separation by the three filters. Photographic 
separation of dye layers produces three separate black and white trans- 
parencies which must be run through the microdensitometer individually. 

Microdensitometry . A scanning microdensitometer was utilized to 

obtain film density measurements of the black and white transparencies. 

This technique permits the rapid scanning density measurement of many 

small adjacent lines in sequence. An Optronics Inc. P-1000 digitizing, 

rotating-drum microdensitometer vas used in this study. This unit scans 

transparencies by way of a rotating drum on which the film is mounted. 

A light source and optical system are used to project a beam through the 

film, and a photo cell mounted below the drum senses the transmitted 

light. The electrical signal representing transmitted light is sampled 

and converted to numerical form on a logarithmic density scale from 0 to 

2. The digital density samples are recorded on magnetic tape with a 

8 

precision of 2 , thus utilizing a 0 to 255 density range compatible with 
the LARS processing system. 

A spatial resolution of 50 um was used, which allowed for a point 
on the film to be resolved which represented about 1.5 cm on the ground. 
Over 300,000 individual density measurements were transformed to a 
numerical data format from a typical 2h by 3 6 mm frame of photography. 



Table 3. Spectral bands from color separation of 
color and color infrared film. 


Separation Filters 



Wratten No. 47b 

Wratten No. 6l 

Wratten No. 23a 

Kodak Ektachrome-X 
and 

Kodak Kodachrome II 

(0.40-0.47wm) 

Blue 

(0.47-0.6lum) 

Green 

(0.59-0.71nm) 

Red 

Kodak Ektachrome 
Infrared 
Type 2443 

(0.47-0.6liim) 

Green 

(0.59-0.71um) 

Red 

(0.68-0.89um) 

Infrared 
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The system has an accuracy of about 2 % in density measurements. 

The recorded film transmission values used are proportional to 
scene reflectance although the exact relationship is not known. Anuta 
and MacDonald (1971) reported that optical film density is a function of 
incoming radiance, film characteristics , filters used, optical system 
characteristics, and the film development process. Densitometry of film 
cannot eliminate characteristics inherently variable in photography, and, 
therefore, is not a calibrated procedure. 

Reformatting and Data Registration . The individual scans were 
stored as separate data sets on tape. In order to analyze the data 
utilizing the information from all three emulsions of the original film, 
the data had to be reformatted in such a manner that each scene point 
was in geometrical coincidence in all three wavelength bands. A soft- 
ware system was developed at LARS (Anuta and MacDonald, 1971) to register 
or "overlay" multiple images of the same scene. The system combined the 
individual scans onto one data set with the image points stored in such 
a way that any scene ground point could be addressed by a line and 
column coordinate pair and a channel pointer. The digital film data 
storage tape was then available for analysis using the LARSYS processing 
system. 

Data points within a field or any rectangular region could be 
identified by extracting line and column coordinates of the corners of 
that area and punching them on computer cards. Each 24 by 36 mm frame of 
photography had about 450 lines and 700 columns of data available for 
analysis. 



Digital Display Imagery . The LARS Digital Video Display Unit was 
used in a preliminary checkout of the data. This device reproduces a 
pictorial representation of the digitized imagery by dividing the data 
up equally into 16 bins for gray level display. A typewriter terminal 
and keyboard are used to control program execution and display functions. 
A light pen is available for selection of coordinate points in the 
imagery. Photocopies of the displayed imagery can be made from a slave 
photocopy unit which is part of the system. 

An initial checkout was made of all 20 frames of digitized photo- 
graphy. It became evident that a problem had been encountered in digi- 
tization of the 1 6 frames of color film since their digital display 
images contained portions which were highly oversaturated. Bright 
portions of the imagery had exceeded the scale of 0 to 255 and had 
"turned over" to low digitization values. The 1 6 frames of digitized 
color film were not used in analysis because of the apparent failure of 
the microdensitometer to accommodate the full range of density values 
represented on the color film. 

Three frames of digitized color infrared film appeared to be of 
high image quality and were chosen for detailed analysis. The light pen 
was used to select a set of data from each frame of digitized photography 
for ground cover analysis . Boundaries were drawn equidistant between 
rows for two or three rows in the frame so that the outlined area could 
be used for an accurate estimation of percent ground cover. Figure 5 
shows the digital display images of the three frames of digitized color 
infrared photography and the areas outlined for detailed study in each 
of them. The imagery for Chalmers plots 9 and 11 is that of the 0.68 to 
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0.89 uni near infrared wavelength band. All three wavelength bands 
(green, red, and near infrared) are shown for Chalmers plot 12. A 
smaller subset area has been outlined in the imagery for Chalmers plot 12. 
This subset area received very detailed study because of the comnlex 
vegetative cover represented therein. Figure 6 shows prints from the 
color and color infrared transparencies taken over these three plots. 

The digitized color infrared photography was obtained from these same 
color infrared transparencies. The high quality of the digitized imagery 
is striking, and one might easily mistake it for the original photo- 
graphic imagery. 

Training Class Analysis - Clustering . Analysis of the digitized 
photography was intended mainly as am attempt to distinguish between the 
components of green vegetation and bare soil in the imagery. This 
involved the determination of separable classes in the data set. The 
statistical pattern recognition analysis system LARSYSAA was used for 
classification purposes. This system uses a maximum likelihood ratio 
based on the Caussian assumption to classify statistically separable 
classes (Fu, Landgrebe, and Phillips, 1969). "Training sample areas" 
must be provided to supply the LARSYSAA processor with statistics for 
recognition of spectral patterns of density values for similar conditions 
in the photograph. A useful technique for selection of training sample 
areas or training classes is the clustering method described by Wacker 
and Landgrebe (1970). The clustering program utilizes a process of 
iteration to minimize the Euclidean distance between relative density 
vectors (in the case of digitized photography) and the "cluster center" 
for a given class. Image points are grouped around these cluster 
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centers such that the overall variance of the resultant "clusters" is 
minimized. 

The clustering program permits calculation of a requested number 
of cluster classes in any desired combination of wavelength bands. For 
the purposes of this study, ten cluster classes were requested from each 
set of data and only the 0.68 to 0.89 urn near infrared wavelength band 
was used for vector analysis. Selection of the near infrared wavelength 
band for clustering purposes was based on the visual observation that 
this band provided better separation of the green vegetation-bare soil 
complex. Ten cluster classes were anticipated to represent sufficiently 
all of the spectrally separable ground features present in a frame of 
digitized photography. 

An area of 100 lines by 100 columns was selected from the center 
of the outlined area on all three frames of digitized photography. This 
area was clustered as described above, and ten spectrally separable 
classes were produced from the imagery from Chalmers plots 9 and 11. 

The clustering program could not distinguish ten spectrally separable 
classes in the imagery from Chalmers plot 12, so it decremented the 
number of classes until it found five cluster centers to which the data 
could be assigned. Checking back to the digital display imagery for the 
three frames of photography, it could be seen that only one cluster 
class from each frame was necessary to characterize the training class 
component for bare soil. Most of the variability in the digitized 
photography was in the green vegetation, and all the remaining cluster 
classes were assigned as training classes for green vegetation. 
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Coordinates of uniform blocks of data representing each cluster 
class were selected by hand from the clustered array map of the cluster- 
ing program. These coordinates were punched on cards, and decks were 
compiled representing training classes for input into the pattern 
classification processor. 

Classification . The smaller subset area outlined on the digital 
display image for Chalmers plot 12 (Figure 5) was chosen for detailed 
analysis and determination of the best wavelength bands for classifi- 
cation. The five cluster classes produced from the clustering of 
Chalmers plot 12 were used for training the pattern classifier. The 
statistics processor of the LARSYSAA pattern classification system was 
used to compute the training statistics for each of the five cluster 
classes. These statistics consist of means and covariance matrices to 
be used by the Gaussian maximum likelihood classifier. 

The coincident spectral plot of means plus and minus one standard 
deviation for the five cluster classes in all three wavelength bands of 
the digitized photography are shown in Figure 7. Cluster classes 1 
through U represent green vegetation while cluster class 5 represents 
bare soil. It can be seen that cluster class 5 (represented by "E" on 
the printout) has a much lower mean in the 0.68 to 0.89 inn wavelength 
band than any of the other cluster classes , and the class representing 
bare soil should not be too hard to differentiate from the classes 
representing green vegetation. The situation in the 0.1*7 to 0.6l pm 
and 0.59 to 0.71 pm wavelength bands is quite different, however. The 
means and standard deviations of cluster class 5 in these two wavelength 
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Figure 7. Coincident spectral plot for 5 cluster classes from digitized photography of Chalmers 
plot 12. 
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cluster classes. Examination of the coincident spectral plot seemed to 
indicate that the 0.68 to 0.89 urn near infrared wavelength band would 
be far better than the other two in discrimination between bare soil and 
green vegetation. 

The LARSYSAA classification processor was then used to classify 
each image point into one of the five defined classes using statistics 
from any desired combination of wavelength bands. Three classifications 
were made using the three wavelength bands by themselves , and a fourth 
classification was made using the statistics from all three wavelength 
bands. Green vegetation classes and display symbols were pooled for 
pictorial displays of com versus bare soil as well as for calculation 
of the percent area occupied by green vegetation. 

Based on preliminary information indicating that the near infrared 
wavelength band was the most useful for classification purposes, the 
other two frames of digitized color infrared photography (Chalmers plots 
9 and 11) were classified using the 0.68 to 0.89 Mm band. 

Results and Discussion 

Relating Estimated Percent Ground Cover to Leaf Area Index 

A regression analysis was performed on the data from the 30 frames 
of photography analyzed for percent ground cover by the point grid sys- 
tem. The data in Table 2 for estimated percent ground cover were plotted 
against leaf area index (LAI). The resulting relationship can be seen in 
Figure 8. It was desired to see if the ground level measurement of LAI 
could be used in this situation to estimate percent ground cover as 
determined from a view above the canopy. If such a relationship exists 
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Figure 8. Percent ground cover plotted against leaf area index for 30 
observations. 
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it would greatly aid the ground based observer whose job it is to collect 
supporting information on ground cover for remote sensing overflights. 

In this problem percent ground cover was considered the dependent 
variable and leaf area index the independent variable. From the graph of 
the relationship (Figure 8) it was felt that a second order term may be 
needed to help explain variations in the dependent variable, so the 
squared term for LAI was also entered into the regression analysis. A 
stepwise regression program was run, and the linear term for LAI was 
entered first into the regression, followed by the squared term for LAI. 
The variation about the mean percent ground cover explained by the re- 
sulting regression was 97.9$« The line of fit is plotted in Figure 8, 
and the regression equation, Y = 2.0972 + 49.3525X - 6.6403X 2 is 

given. For this particular situation, leaf area index appears to be 
useful in explaining variations in percent ground cover. 

It is evident from Figure 8 that as LAI increases beyond a certain 
point (LAI * 3 in this case) percent ground cover does not increase very 
greatly. A ground cover of 100$ was never observed from the low alti- 
tude photography. Cuellar (1971) found it difficult to estimate cover 
from the ground as ground cover approached 100$. If a prediction equa- 
tion like the one calculated were available for different row widths and 
plant populations of corn, estimation of cover from LAI measurements 
would not be unreasonable, even as cover approached 100$. 

Results of Classifications Using Three Color Separations 

Display maps of the classified subset area for the digitized photo- 
graphy from Chalmers plot 12 are shown in Figure 9. Four display maps 
are shown, representing classifications in the green, red, and near 




46 


infrared wavelength bands separately and all three bands in combination. 
Green vegetation is displayed as a "C" and bare soil is displayed as a 

The classification using the near infrared wavelength band was 
the best one for determination of percent ground cover in the subset 
area of Chalmers plot 12. The classification display map of this area 
in the near infrared wavelength band represents most closely the actual 
ground cover situation when compared with the original photography 
(Figure 6) and the digital display imagery (Figure 5). 

The tabulated results of percent ground cover for the subset area 
of Chalmers plot 12 as determined by the four classifications are listed 
in Table 4. It earn be 3een that percent ground cover was overestimated 
in the classification using the red wavelength separation, while it was 
slightly underestimated in the classification using the green wavelength 
separation. The classification using all three wavelength bands in 
combination very greatly underestimated percent ground cover as is evi- 
dent from the display map in Figure 9* as well as the results in Table 4. 
About one-third of the points in the map display for the multi-band 
classification were thresholded and are displayed as blanks. A threshold 
level of 0.1 % was applied to prevent assignment of a point to a par- 
ticular class when membership in that class was not highly probable. 

Classifications using the near infrared wavelength band were made 
from the digitized photography for Chalmers plots 9 and 11. A larger 
area from Chalmers plot 12 was also classified for percent ground cover 
determination. Figure 10 shows the display maps for these three plots. 
From the display map for Chalmers plot 11 it can be seen that even some 
green weeds in the plot were classified as green vegetation. The 


Table 4. Percen 
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classification results from these three plots and the corresponding 
leaf area indices are given in Table 5. 

Evaluation of Color and Color Infrared Film 
in Determining Ground Cover 

The color and color infrared photographs (Figure 6) show the 
differing usefulness of each in estimating ground cover on a dark soil. 
The color photographs are not as valuable as the color infrared photo- 
graphs mainly because of the difficulty in distinguishing leaves in the 
shaded areas between the rows. This is in agreement with the findings 
of Hoffer ( 1967 ) who noted that increasing canopy on dark soils resulted 
in higher response on infrared film. Gerbermann, Gausman, and Wiegand 
(1969) likewise, found that the influence of crop shadows and sunlit 
furrows could be detected on color infrared film. 

On the light colored Russell plots color film seemed to be the 
most useful for estimating ground cover when the percent cover was less 
than about 75 % and the soil was in a typical dry surface condition. If 
the surface soil of the Russell plots was very moist and the ground 
cover was greater than about 75 %, however, color infrared film seemed to 
be best for estimating ground cover. Often when a corn canopy reaches 
75 % ground cover, the evaporation from the surface soil is lessened, and 
the surface soil appears darker. If this is the case, color infrared 
film would be preferable in either a dark or light soil background 
situation in which the percent ground cover is greater than about 75 %- 
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Heed for More Quantifiable Information 

Differences in exposure, development, and manufacturing of photo- 
graphic films prevent simple quantitative calibration of photography to 
reflectance values. Although reflectance calibration of photographic 
film has been attempted by the use of gray scale step wedges and ground 
reflectance panels (Silvestro, 1969), it is not at all practical for the 
type of experiment described here. Very little information as to the 
actual reflectance values could be obtained from the photography taken 
in this study. 

Microdensitometry and digitization of photographic film can be a 
precise and quantitative process with known allowable error. Unfortu- 
nately, however, the photographic film from which the data is obtained 
is subject to many sources of variation as listed above. Digitized 
photography is only as good as the original photograph. 

Photography is limited to the wavelength region of about O.U to 
0.9 um. There may very well be other wavelength regions which may be more 
useful in estimating ground cover and which may be determined only by 
non -photographic sensors. Multispectral airborne scanners and portable 
field spectro radiometers may be able to cover the entire reflective wave- 
length range while obtaining more quantifiable information on ground 
cover. 


Summary and Conclusions 

Color infrared film is potentially useful for ground cover deter- 
mination from low altitude photography. The near infrared dye layer of 
this film was the most useful for separation of the green vegetation and 
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bare soil components of the com canopy. Digitization of low altitude 
photography of plant canopies provides a valuable method for computer 
assisted determination of percent ground cover. 

Color film is of value in ground cover determination from low 
altitude photography of com canopies if the soil background is light 
and the ground cover is less than about 75 %. For dark soil backgrounds 
and ground covers of greater than 75 % in com canopies, color infrared 
film would be more valuable. 

Percent ground cover as determined from low altitude photography 
can be related to leaf area index of a com canopy for a given row width 
and plant population. Leaf area index is easy to measure from the 
ground and gives more information as to the actual crop biomass than 
does percent ground cover. 

Low altitude photography of corn canopies is useful for spatial 
characterization of ground cover, but quantifiable spectral information 
is limited due to variability between frames of photography. Micro- 
densitometry and digitization of color and color infrared film provide 
quantitative techniques for analyzing density differences vhich may be 
related to components of green vegetation and bare soil in the scene. 

Since the near infrared color separation of color infrared film 
seemed to be the most useful in ground cover determination, it may be 
Just as valuable to use black and white infrared film filtered to the 
near infrared for photography of com canopies. Color infrared film is 
still of value even without having its three dye layers photographically 
separated. 
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CHAPTER IV 

AIRBORNE MULTISPECTRAL SCANNER DATA FROM CORN CANOPIES 
Problem Definition 

Airborne multi spectral scanner data provide an extended view of the 
reflective properties of plant canopies. Information may be obtained in 
the wavelength range from O.U urn to 2.6 um in this manner. Comparisons 
of scanner data values in the various discrete wavelength bands or 
channels of the multispectral scanner allow for investigation of the 
spectral differences between various plant canopy covers. 

Kristof and Baumgardner (1970) found a ratio technique to be useful 
in relating multispectral scanner data values to extent of ground cover. 
In particular, the ratios of scanner data values from the near infrared 
and chlorophyll absorption regions of the spectrum appear to be of great 
value in expressing changes in ground cover. These ratios can be related 
to a ground based measurement of ground cover such as leaf area index. 
Prediction equations can then be formulated for estimation of leaf area 
index from ratios of multispectral scanner data values. 



Procedures 


Multispectral Scanner Data Collection 

Two flight missions by the University of Michigan's airborne multi- 
spectral scanner were used in this study. The flight path for both 
missions was a north-south flightline directly over the Purdue University 
Agronomy Farm's Weather Station. Corn plots of various ground covers 
planted on a Russell silt loam soil were located in the flight path in 
close proximity to the weather station. The weather station was the 
location of placement of a set of eight ground reflectance panels whose 
use in standardizing scanner response over time will be described later. 
Figure 11 shows the placement of the panels with respect to the weather 
station. 

July 12 and July 21, 1971 were the dates for the two scanner over- 
flights. The July 12 flight was at 1401 hours EST while the July 21 
flight was at 1305 hours EST. Both overflights were at an altitude of 
305 m and both took place on clear, cloud-free days. 

The University of Michigan's mult ispect ral scanner collected spectral 
data in the 0.46 to 2.60 pm reflective wavelength range, as well as the 
9.3 to 11.7 pm thermal infrared wavelength range. Only the reflective 
wavelength bands or channels were used in this study. Table 6 lists the 
12 channels of data collected by the University of Michigan scanner and 
the corresponding wavelength bands. 

The aircraft scanner data were recorded in analog form and converted 
to a digital format for analysis by the LARSYS processing system. Two 
flightlines of interest were digitized at a sample rate of 450 samples 
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Figure 11 . Ground reflectance panel placement at the 
Purdue University Agronomy Farm Weather 
Station, August 17 , 1971 . 
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Table 6. 

The 12 channels and corresponding wavelength bands for the 
University of Michigan multi spectral scanner. 

Limits of Spectral Bands (pm) 

Channel 

Lower 

Upper 

Wavelength Region 

1 

0.46 

0.49 

Visible 

2 

0.48 

0.51 

Visible 

3 

0.50 

0.54 

Visible 

4 

0.52 

0.57 

Visible 

5 

0.54 

0.60 

Visible 

6 

0.58 

0.65 

Visible 

7 

0.61 

0.70 

Visible 

8 

0.72 

0.92 

Near Infrared 

9 

1.00 

1.40 

Near Infrared 

10 

1.50 

1.80 

Near Infrared 

11 

2.00 

2.60 

Near Infrared 

12 

9.30 

11.70 

Thermal Infrared 
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per scan line of data. This meant that every line of analog scanner 
data was sampled in the digitization process. Each resolution unit 
represented an area about 1 m by 1 m on the ground. The digitized data 
were reformatted into a form compatible with data processing programs at 
LAFS. 


Digital Display Imagery 

The LARS Digital Video Display Unit was used to display the digi- 
tized scanner data imagery in the twelve wavelength bands of the multi- 
spectral scanner. The digital display images for the July 12 scanner 
flight in the twelve wavelength channels for the area near the Purdue 
University Agronomy Farm's Weather Station are shown in Figure 12. These 
display images were obtained from black and white Polaroid photographs 
taken from the photocopy unit which is part of the digital display sys- 
tem. The images are composed of sixteen gray levels which were produced 
by dividing all the data points in the displayed area into sixteen equal 
bins of data. Even though the images look like normal panchromatic 
photographs , it must be remembered that each image represents the rela- 
tive response values only in the wavelength band corresponding to that 
channel (Figure 6). 

The ground reflectance panels can be seen in the center of the 
imagery. These panels are 6.1 m by 12.2 m standard reflectance canvases 
painted with special reflecting paints. The gray reflectance panels have 
a relatively constant percent reflectance throughout the reflective wave- 
length range, and are useful for visual comparisons with ground features 
in the nearby area. The reflectance panels for the July 12 flight were 
arranged in the same order as in Figure 11 (taken August 17) except that 



Channel 1 (0. 46-0. 49pm) Channel 2 (0 . 48-0 . 51ym) Channel 3 (0. 50-0. 54pm) 



Channel 4 (0 . 52-0 . 57pm) Channel 5 (0 . 54-0 . 60pm) Channel 6 (0 . 58-0 . 65pm) 



Channel 7 (0 . 61-0 . 70pm) Channel 8 (0 . 72-0. 92pm) Giannel 9 (1.0-1. 4pm) 



Channel 10 (1.5-1.8um) Channel 11 (2. 0-2. 6pm) Channel 12 (9 . 3-11. 7ptn) 


Figure 12. Digital display imagery of Russell plots and panels^ July 12 
scanner flight. 
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the three color panels for July 12 were arranged green-blue-red from east 
to vest. 


Locating Plots and Obtaining Mean Scanner Data Values 

The approximate area of the Russell plots was outlined on the digi- 
tal display screen using the light pen, and the coordinates of the area 
were recorded by the field selection processor. This area was then 
clustered utilizing the data from all twelve wavelength bands of the 
multispectral scanner with ten cluster classes being requested. The 
resulting cluster array map was of value in locating the twelve sub-plots 
within the Russell plots. Sub-plots that were planted on the same date 
and were of approximately the same ground cover were grouped into the 
same cluster classes by the clustering program, thus facilitating the 
selection of data blocks representing each sub-plot. 

Each sub-plot was slightly larger than U by U ground resolution units 
or data points. A square block of four data points was selected from the 
northeast portion of each sub-plot in order to eliminate any shadowing 
effects from bordering sub-plots. Statistics were then obtained from 
these data blocks for all of the sub-plots on both dates of scanner over- 
flights. 

At the time of the July 21 scanner overflight, photographic data 
were being obtained from the LARS Hi-Ranger truck over Russell plot 2. 

The location of the Hi-Ranger truck in relation to the Russell plots can 
be seen in Figure 13. This print is from a black and white photograph 
taken by the scanner airplane, coincident with the scanner data. It 
became evident upon cluster analysis of the Russell plots for the July 21 
scanner overflight that the Hi-Ranger truck was having an effect on the 


(Overlays page 60) 






Figure 13. Black and white scanner coincident photograph filtered to the red wavelength region, 
taken July 21, 1971 over the Purdue University Agronomy Farm Weather Station. 
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spectral response of three of the sub-plots. For this reason, Russell 
plots 2, 3, and U were discarded from further analysis of the scanner 
data from July 21. 

The LARSYS statistics processor was used to obtain mean scanner data 
values for the sub-plots for the two flight dates. These mean values 
were obtained without any attempt made to calibrate scanner response 
between flights. The normal procedure used at LARS for obtaining statis- 
tics from multispectral scanner data was used, which takes into account 
zero reference drift or the drift in the scanner system response to a 
black or zero signal level. More about calibration procedures for multi- 
spectral scanner data will be presented in the next chapter. The mean 
scanner data values obtained for the twelve sub-plots in the eleven 
reflective wavelength bands for the two dates are given in Table 7. 

Ratio Techniques for Relating Scanner 
Response Ratios to Leaf Area Index 

JCristof and Baumgardner (1971) used the ratio of scanner data values 
in the 0.58 to 0.62 pm wavelength band to the relative reflectance in the 
0.80 to 1.00 pm wavelength band to observe increases in ground cover. 
Considerations of the reflectance properties of green vegetation and soil 
would seem to favor a ratio of these two wavelength regions for character- 
izing ground cover. Chlorophyll absorbs at 0.65 pm in the visible wave- 
lengths , and it is reasonable to assume that as a crop canopy becomes 
more dense it would absorb more strongly (and hence, reflect less) in this 
region. The 0.72 to 1.3 Pm wavelength range is the region in which little 
absorption takes place in green vegetation and in which almost all in- 
cident light is transmitted or reflected. Multiple leaf layer experiments 


Table J. Uncalibrated mean scanner data values for 12 Russell plot 3 in 11 reflective wavelength 
bands for two flight dates. 
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have shown that reflectance in this 0.72 to 1.3 um wavelength region 
increases with increasing leaf layers (%ers et al. , 1966). Therefore, 
the reflectance of a crop canopy may be expected to increase in the 0.72 
to 1.3 ym wavelength range as green vegetative cover increases. Soil 
reflectance generally increases in going from the red wavelength region 
to the near infrared wavelength region, the amount of increase for a 
particular soil depending mainly on the soil moisture content. 

The scanner wavelength band configuration was changed from that at 
the time of the study of Kristof and Baumgardner (1971) » the chlorophyll 
absorbtion band now falling in the 0.6l to 0.70 ym wavelength band of 
channel 7. Two near Infrared wavelength bands were now available in the 
region of high reflectance of green vegetation; the 0.72 to 0.92 pm 
wavelength band of channel 8; and the 1.0 to 1.4 ym band of channel 9. 

Ratios of scanner data values (Table 7) were calculated by dividing 
the scanner data values from channels 8 or 9 by the scanner data values 
from channel 7 (hereafter referred to as the ratios of channels 8/7 or 
9/7). These ratios were then plotted against leaf area index, as measured 
at the time of the overflights. The ratios of channels 8/7 and 9/7 for 
the two flight dates and the corresponding LAI's for all of the sub-plots 
are listed in Table 8. Stepwise multiple regression analysis was run on 
the data to see if either of the two ratios of near infrared to red wave- 
length scanner data response could be related to leaf area index. 
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Table 8. Ratios of uncalibrated scanner data values and leaf area 
indices for 12 Russell plots for two flight dates. 

July 12, 1971 Multispectral Scanner Flight 

Russell Ratio Ratio 


Plot No. 

LAI 

Channels 8/7 

Channels 

1 

0.01 

0.56 

0.76 

2 

2.89 

1.63 

1.97 

3 

3.48 

1.83 

2.18 

1+ 

1.30 

1.1k 

1.U9 

5 

0.01 

0.58 

0.80 

6 

2.34 

1.51 

1.78 

7 

1.11 

1.20 

1.52 

8 

0.04 

0.63 

0.9!+ 

9 

3.10 

1.68 

2.02 

10 

1.17 

1.23 

1.59 

11 

2.93 

1.62 

1.96 

12 

3.82 

1.75 

2.09 


July 21, 1971 Multispectral Scanner Flight 

0.9 1 * 

Missing 
Data 

0.94 
2.08 
1.88 
0.98 
1.95 
1.85 
2.08 
2.08 


1 

0.11 

1.21 

2 

3.00 

f 

3 

U.28 

Missing 
| Data 

U 

2.31 

l 

5 

0.12 

1.09 

6 

3.57 

2.71 

7 

2.15 

2.39 

8 

0.17 

1.15 

9 

3.38 

2-58 

10 

2.10 

2.42 

11 

3. 45 

2.66 

12 

U .06 

2.67 
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Results and Discussion 

Visual Comparisons of Digital Display Imagery 

The photographs obtained from the LARS Digital Video Display Unit in 
the twelve wavelength bands of the multispectral scanner are shown in 
Figure 12. The Russell sub-plots are rather small to be able to pick out 
each one individually, but the low ground cover plots do stand out and 
serve as reference points. The speckled appearance of the imagery in 
channels 4 and 7 was caused by problems in digitization of the multi- 
spectral scanner data, but this problem does not seem to have affected 
any of the ground features of interest. The bulk corn plots to the east, 
west, and south of the reflectance panels are of special interest when 
viewed in all wavelength bands. These bulk corn plots were planted near 
May 1, and had a very dense cover by this date. 

In channels 1-7 the bulk com fields have a very low reflectance 
when compared to the gray ref lectern ce panels. It can be seen that the 
bulk com fields are reflecting even less than the black reflectance panel 
except in the peak of the green wavelengths in channel 4 ( 0.52 to 0.57 um). 
This black reflectsuice panel has a mesisured percent reflectance of not 
more than h% throughout the reflective wavelength range according to DK-2 
spectrophotometer measurements. If the black reflectance panel is indeed 
reflecting less than h% in the field, then the com canopies are appar- 
ently reflecting even less than 4$ in all but channel 4 of the visible 
wavelengths. This would tend to be in agreement with Knipling (1969) who 
estimated that the visible wavelength reflectance from a complete canopy 
may be 3 to 5 percent reflectance. 
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In channels 8 and 9 , the bulk corn fields are reflecting somewhere 
between the reflectance of the white and the light gray panels. DFC-2 
spectrophotometer measurements of the white and light gray panels indicate 
that they reflect 56 % and 27 % respectively in the wavelength region of 
channels 8 and 9. If these panel reflectances are similar in the field, 
the bulk corn canopies are reflecting between 27 % and 36 %. This is also 
in agreement with Knipling's estimate that the near infrared reflectance 
from a complete canopy would be about 35 %. 

The reflectance of the bulk com plots decreases in channel 10 from 
the high reflectances in channels 8 and 9, and then decreases in channel 
11 from the reflectance in channel 10. The reflectance of the corn plots 
follows the same general response as that of green vegetation studied in 
the laboratory (Hoffer and Johannsen, 1969 ). The reflectance panels do 
hot have standard thermal properties , but it can be seen from the imagery 
in channel 12 that all of the panels seem to be quite warmer than their 
surroundings, especially the relatively cool cornfields. 

Scanner Response Ratios 
and Leaf Area Index for Two Flight Dates 

The ratios of scanner data values in channels 9/7 for the July 12 
and July 21 flight dates were plotted against leaf area index (Figure lU); 
The points seemed to indicate that a rather good relationship existed 
between the ratio of scanner data in channels 9/7 and LAI. A stepwise 

n| 

multiple regression analysis indicated that 96 . 8 % of the variation about 

the mean LAI could be explained by the second order regression equation, 

* 

Y 3 1.0751 - 2.8756X + 1.9M2X 2 , where X is equal to the ratio of scanner 

data values in channels 9/7. This indicated that ratios of multispectral 




0 s July 12 Flight 
X= July 21 Flight 


L r*l.07Si-2.8756X+ I.9442X' 

^ . / Channel 9 \ 

whew X» Ratio (•— ■ — 

\Chonnel 7 / 


R 2 - 0.968 


' Scanner data value in channel 9 


Scanner data value in channel 7 

Leaf area index versus the ratio of scanner data values in 
channels 9/7 for two flight dates. 
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scanner data may, as hoped for, "be related to sane ground based measure- 
ment of ground cover. It was interesting to note that the ratios of 
scanner data values for both dates compared quite well even without any 
calibration of scanner response between the two dates. 

It was expected that the ratios of scanner data values in channels 
8/7 should also show a good relationship with LAI. The ratios of scanner 
data values in channels 8/7 for the July 12 and July 21 flight dates were 
plotted in Figure 15. It very soon became obvious that two separate 
relationships existed for the two multispectral scanner flights. Step- 
wise multiple regression analysis indicated that 98.1$ of the variation 
in LAI could be explained by the second order regression equation, 

Y * -0.^875^ + 1.26899X 2 , where X equals the ratio of scanner data values 
in channels 8/7 for the July 12 scanner flight. For the July 21 scanner 
flight a linear relationship was indicated which was able to account for 

-93.0$ of the variation in LAI. The regression equation for July 21 was 

A 

Y = -2.U611 + 2.182U5X, where X equals the ratio of scanner data values 
in channels 8/7 for the July 21 scanner flight. 

A relationship did seem to exist between the ratio of scanner data 
values in channels 8/7 and LAI , but the same relationship did not hold 
between flight dates. Apparently a change had taken place between 
July 12 and July 21 in the way in which the multispectral scanner was 
recording data in channel 8. This indicated that raw scanner data could 
be highly unreliable in any attempts to use a ratio technique over time 
(between flight dates ) . 



0 = July 12 Flight 
X» July 21 Flight 


Y* -0.48754 +■ I.26899X- 


where X s Ratio 


r Channel 8 
^Channel 7 


R 2 = 0.981 


/,Y=-2.46III + 2.I8245X 

u v/ « * /Channel 8 \ 
Wh6re X=Ro, '° \ Channel 7 ) 

R 2 * 0.930 


Scanner data value in channel 8 


Scanner data value in channel 7 

Figure 15. Leaf area index versus the ratio of scanner data values in 
channels 8/7 for two flight dates. 
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Problems in Evaluating Scanner Data Over Time 

Comparisons of multispectral scanner data between flight dates have 
always been difficult because of the many variable factors involved. 
Weather and atmospheric conditions , scene illumination intensity sis a 
function of wavelength, time of day, and angle of illumination can always 
be counted on to complicate comparisons between multispectral scanner 
flights . 

Not only are there natural phenomena to contend with, but there are 
also many involving the scanner system itself. Data values for the same 
ground target have been observed to change from one side of the flightline 
to the other and from the beginning of a flightline to the end. Changes 
in scanner response over time within the same flightline may occur due to 
drift in zero level reference as well as gain changes in the system. Gain 
changes are often made in one channel and not in another, thus it becomes 
difficult to make any comparisons between channels over time. 

Many potential applications of remote sensing depend on the 
ability to view repeatedly a target of interest and characterize the 
spectral properties of that target over time. Ground cover sensing 
is certainly an area in which this ability is needed. A more 
sophisticated calibrated airborne spectrophotometer is needed to 
increase the capability to compare scanner data over time. 
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Summary and Conclusions 

Multispectral scanner data from corn canopies provide more quanti- 
tative spectral information on ground cover than does photography. Ratios 
of uncalibrated scanner data can be used within a given scanner mission 
to estimate changes in vegetative cover as measured by leaf area index. 

The peak reflectance wavelengths for green vegetation (0.72 to 1.3 ym) 
and the region of strong chlorophyll absorption (0.65 ym) seem to provide 
the best information on changing ground cover. 

Regression equations can be evolved relating leaf area index in 
com canopies to the ratios of scanner data values from channels 8 and 
9 (0.72 to 0.92 pm and 1.0 to 1.4 ym) to scanner data values from 
channel 7 (0.6l to 0.70 um). This can only be done within a given 
flight mission unless the investigator is certain that there were no 
changes in scanner system gain between flights. 

Visual observation of digital display imagery of the ground reflec- 
tance panels and nearby cornfields indicated that dense com canopies 
reflected less them 4# in. most visible wavelength bands. In the 0.72 to 
1.4 ym wavelength range, the dense com canopies reflected between 27# 
emd 56 # when compared to the ground reflectemce panels. These values 
were within the expected ranges of reflectance estimated by Knipling 
(1969). 
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CHAPTER V 

NORMALIZED SPECTRAL RESPONSE OF CORN CANOPIES FROM 
AIRBORNE MULTISPECTRAL SCANNER DATA 

Problem Definition 

Airborne multispectral scanner data allow for examination of the 
spectral differences between various plant canopy covers. Ratios of 
scanner data response can be related to the ground based measurement of 
leaf area index. It is desirable to be able to compare results from 
more than one flight date. In this way the theorized relationships 
between ratios of scanner data values and leaf area index can be tested. 

Variations in scanner system response between flight dates have 
prevented direct comparison of scanner data over time. Internal cali- 
bration standards within the multispectral scanner and reference to 
ground reflectance panels permit normalization of scanner response 
between flight dates. The objective of this study was to attempt to 
normalize the spectral response of corn canopies as detected by an air- 
borne multispectral scanner. 
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P rocedures 

Internal Calibration - Black Level and Calibration Lamp 
The University of Michigan multi spectral scanner records three 
calibration sources for every scan line of data in the reflective wave- 
length region. These three calibration values are stored on the 
Laboratory for Applications of Remote Sensing (LARS) aircraft data 
storage tape and may be called upon by the data analyst to effect inter- 
nal calibration of scanner response. The three calibration sources are: 

1. CO, a black or dark level reference which gives a measure of black or 
zero signal level, 2. Cl, a constant light source filtered to approximate 
closely the spectral energy curve of solar illumination, and 3. C2, a 
sun sensor which measures the solar illumination incident upon the top of 
the aircraft. These three sources have the potential to permit compari- 
son of multispectral scanner data collected at different times. The 

j 

black level source, CO, and the precision calibration lamp source, Cl, 
should remain constant with the passage of time, but the sun sensor 
source, C2, will vary according to solar illumination. 

The black level, CO, and precision calibration lamp, Cl, are 
completely internal standards whereas the sun sensor, C2, depends on 
conditions external to the scanner. The sun sensor records the solar 
illumination at aircraft altitudes, which may or may not be the same as 
the solar illumination of the target. Sun sensor calibration becomes 
particularly undependable on cloudy days when clouds may be shadowing 
the sun sensor but not the ground target. For purposes of this study, 
only the CO and Cl calibration sources were used. 
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The calibration lamp values recorded with the scanner data repre- 
sent a precise luminous energy source which varies only with the lamp 
current, lamp age, and lamp filter used. Use of the calibration lamp 
values compensates for gain changes in the scanner system. 

The conversion of analog scanner data to digital form is performed 
at LARS with the aim of utilizing the maximum dynamic range of the data 
system. The amplifier gain changes necessary to utilize this dynamic 
range and the characteristics of the scanner system itself, result in the 
fact that the data values in each channel are unrelated to those in other 
channels. Within a channel, however, data values are directly related 
to scene reflectance by some unknown function. 

A linear calibration function is used in the LARSYS processing 
system (a computer software system developed at LARS for analysis of 
aircraft scanner data) which delivers data to the statistics processor. 
Phillips (1969) describes the implementation of this function in the 
calibration procedure. The calibration function is: 

°D * A C ‘ D T + B C 
where 

Dp * Data delivered to the statistics processor 
D t = Data stored on the aircraft data storage tape 
A c and B c = Calibration values 

A c and are calculated for each scan line of data and are used through- 
out the entire scan line for which they have been calculated. For this 
study, Ap and B^ were set to the values for the black level and the 
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calibration lamp, thus utilizing a two point calibration capable of 
removing lew frequency bias level drift and amplification changes from 
the system. A c and B c could have been set to any of seven combinations 
of the three calibration sources, but this particular combination was 
chosen because the system drift rad gain were known to vary from one 
end of the flightline to another. 

The calibration values for A c and B^ in the above equation were 
extracted from the values for CO and Cl in the identification record for 
the particular flightline. The average CO and Cl values for channels 1 
through 11 were calculated and recorded for the July 12 and July 21 , 1971 
flight dates (Table 9). These CO and Cl values were used to obtain two 
point internal calibration of the scanner data response from the Russell 
silt loam test plots and from the set of ground reflectance panels. 

External Calibration - Ground Reflectance Panels 

Locating Panels and Obtaining Mean Scanner Data Values . Eight 
ground reflectance panels were available on the two flight dates 
(Figure 11). The set consisted of 5 gray level panels of nominal re- 
flectances of 8$, 1 6 %, 32% * and 6k% as well as red, green, and blue 
panels. These panels served as a form of external calibration providing 
a ground to aircraft link capable of removing the effect of atmospheric 
scattering (Silvestro, 1969)* Use of the panels allows estimation of 
gain correction factors for approximation of actual scene reflectance in 
each channel. Hasell and Larsen (1968) describe the use of these eight 
reflectance panels in calibrating the output of the University of Michigan 
multi spectral scanner. 
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Table 9. 

Internal calibration values for two 
July 12, 1971 Run No. 71029401 

flight dates. 

Channel 


CO Value 

Cl Value 

1 


24.36 

71.08 

2 


28.27 

105.90 

3 


23.25 

94.83 

4 


29.86 

66.6 3 

5 


22.05 

76.64 

6 


21.64 

92.50 

7 


28.37 

90.10 

8 


19.64 

79.14 

9 


23.92 

33.31 

10 


19.50 

68.89 

11 


23.84 

93.54 


July 

21, 1971 Run No. 71054401 


Channel 


CO Value 

Cl Value 

1 


21.87 

55.60 

2 


21.79 

78.51 

3 


23.55 

48.40 

4 


18.24 

63.75 

5 


26.88 

66. 6o 

6 


21.79 

75.04 

7 


28.72 

75.40 

8 


27.11 

68.63 

9 


23.51 

26.10 

1° 


20.36 

33.49 

11 


22,47 

40.53 
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Calibration to ground reflectance panels permits normalization of 
scanner data to scene reflectance when the area of interest is in 
environmental proximity to the reflectance panels. Environment al proxi- 
mity in this case means an area of the same illumination, the same sun 
angle, the same aircraft altitude, and the same atmospheric conditions 
as the area from which scanner data are collected for ground reflectance 
panels . 

The approximate locations of the reflectance panels in both flight- 
lines were determined by outlining the areas with the light pen on the 
digital display screen. These areas were then clustered using the data 
from channels 2,4,6,8,10, and 12. Ten cluster classes were requested 
in each area, but the clustering program decremented the number of cluster 
classes until only the same number of cluster classes as panels remained. 
The coordinates of the ground reflectance panels were then obtained from 
the cluster array map. Each panel was represented by about 100 data 
points. 

The LARSYS statistics processor was used to obtain internally 
calibrated mean scanner data values for the panels for both flight dates. 
The values for CO and Cl in the identification record for each run were 
used for two point calibration of the scanner response in order to reduce 
variations caused by gain changes and drift in scanner response through- 
out the flightline. For the purposes of this study, all CO values were 
set equal to zero, and Cl values were calculated by subtracting the iden- 
tification record values of CO from the identification record values of 
Cl. This was done in order to set the black level or zero radiation 
standard to zero response, since the black level standard always has the 
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lowest response of any object viewed by the scanner. The dynamic range 
between the identification record values of CO and Cl was thus retained, 
the only difference being the scaling-down of the range so that black 
level response was set to a scanner response of zero. From the resulting 
scanner data values for the eight ground reflectance panels for the two 
flight dates , it can be seen that some of the values for the lighter 
panels saturated or "turned over" because of failure of the scanner system 
to respond correctly to these high reflectance values (Table 10 ). The 
saturated data values could not be used in analysis. 

The scanner data values for the ground reflectance panels (Table 10 ) 
were later used in relating scanner response to actual scene reflectance. 
The next step in the procedure involved the laboratory measurement of the 
percent reflectance of the ground reflectance panels. 

Laboratory Measurements of Spectral Reflectance of Panels . It is 
assumed that the ground reflectance panels behave as perfectly diffuse or 
Lambertian reflectors of incident illumination; that is, they exhibit a 
uniform spatial distribution of radiance, independent of the geometry of 
illumination. Another assumption which had to be made was that laboratory 
DK-2 spectroreflectometer measurements of percent reflectance could be 
related to percent reflectance in a field situation. No field spectro- 
radiometer was available to measure actual directional reflectance of the 
panels in the field so the DK-2 spectroreflectometer was used to char- 
acterize the reflectance of the eight panels. 

Differences exist between the DK-2 spectroreflectometer and field 
or airborne spectroradiometers in the way in which they measure percent 
reflectance. In the DK-2 spectroreflectometer, illumination is normal to 


Table 10. Scanner data values for ground reflectance panels, calibrated to CO and Cl 


80 


CO 

Os 

CO 


O 

CM 

H 

H 

t- 


OJ 

CO 

f — 1 

LA 

CM 

-it 

OJ 

^t 


VO 

• 

• 

* 

• 

• 

• 

• 

• 

« 

• 

CM 

VO 

On 

H 

On 

la 

On 

CM 

m 

On 

CM 

CM 

O 


ITN 

vo 

VO 

ITN 

CO 

CO 

CM 

CM 

CM 






rH 

CM 




O 

LA 

LA 

H 

VO 

ON 

O 

O 

.it 

c 

OJ 

H 


-it 

IA 

CO 

vo 

CM 

rH 

vo 

0 

0 ) 

• 

• 

# 

• 

• 

• 

• 

• 

• 

• 

^t 

O 


-it 

ON 

t- 

On 

-it 

H 

On 

u 

LA 

CM 

OJ 

-it 

H 

CO 

vo 

-it 

CO 

CO 

0 

•H 

CM 

CM 

rH 

H 




H 

CM 


H 

O 

-it 

ON 

CM 

o 



CO 

00 

CO 

On 

On 

1A 

CM 

-it 

-it 

CO 

Tfl 

CO 

OJ 

CO 

O 

LA 

IA 

CO 

VO 

OJ 

CO 

<1)1 

• 

• 

' • 

• 

• 

• 

• 

• 

• 

• 

ft; 

0 

LA 

t- 

On 

O 

On 

00 

t— 

CO 

H 


vc 

co 


CO 

vo 

vo 

rH 

CO 

t- 

00 







«H 

CM 

rH 

H 

CM 


«S” 

-it 

Ol 

On 

c- 

vo 

CM 

OJ 

O 

CM 


CO 

LA 


H 

-it 

LA 


t- 


O 

• 

• 

• 

• 

• 

• 

• 

• 

• 


LA 

-=t 

CO 


t- 

CM 

O 

-it 

CO 

LA 

CO 

t- 

co 


LA 

IA 

CM 

. -it 


CO 

rH 



>> 

-it 

CO 

-it 

rH 


CM 

CO 

CO 

O 

ON 

#* 

a? 


IA 

VO 

O 

-d* 

O 

VO 

vo 

On 

h- 

p 

P 

• 

• 

* 

• 

• 

• 

• 

• 

• 

- • 


O 

O 

8 

iH 

CM 

LA 

-It 

H 

LA 

vo 

rH 

rH 

to 

•H 

& 

Vt 

co 

On 

H 

H 

LA 

vo 


vo 

CM . 

-it . 

h- 


s 












§ 

aS 

f-4 

8 

CO 

vo 

LA 

& 

LA 

CO 

OO 

rH 

CM 

8 

O 

On 

VO 

rH 

0 

O 

• 

• 

. • 

• 

• 

• 

• 

« 

* 

• 

CO 

vo 

vo 

CO 

t- 

CM 

CO 

s 

rH 

vo 

00 



LA 

CM 

0 

On 

H 

CM 

-It 

LA 

ON 

H 

VO 

H 

CM 

H 


H 

rH 

H 





H 












OJ 

H 


g 1 

O 

V* 

CVJ 

oo 


0) 


CM 

-P 

LA 

CO 

VO 

rH 

On 

-it 

On 

CO 

VO 

d 

H 

CM 

H 

On 

IA 

CO 

-it 

rH 

• 

U 

• 

• 

• 

• 

« 

• 

• 

• 

CO 

P 

O 

H 

CO 

CO 

CO 

vc 

CM 

CM 

CM 

P 

CO 

t- 

0 

H 

0 

VO 

CO 

-it 

CM 

<tf 

CM 

H 

H 

CM 

CM 



H 


(0 



T* 










*1 

0) 

0) 

<u 

0) 







P 

p 

p 

p 

vo 

rH 

rH 

-it 

H 

CO 

P 

«5 


d 

d 

CM 

t- 

CM 

O 

vo 

CM 

O 

b 

H 

b 

b 

• 

* 

• 

• 

« 

• 

p 


p 

3 

rt 

H 

On 

t- 

H 

ON 


p 

p 

p 

p 

CO 

CM 

H 

LA 

CO 

-it 

a 

oJ 

«$ 

a) 

fvj 

CM 

CM 

H 

rH 

CM 

vo 

to 

co 

CO 

co 








OJ 


on 


la 


no 


co 


o 

i— J 


23U.W* 1U8.01 102.U3 72.58 71.18 225.72 227.72 230.76 



Table 10, cont 


81 



CO 

-st 

CM 

h- 

CO 

On 

CM 

Os 

o 

h- 

t- 

0) 

CVI 

CO 

t- 

o\ 

CO 

H 

t- 

CO 

CO 

vo 

co 

JS 

• 

• 

• • 

• 

• 

• 

• 

• 

. • 

• 

• 

H 

-st 

H 

on 

On 

t a 

H 1 

LA 

on 

CM 

vc 

CO 


CM 

CM 

o 

H 

vc 

h- 

t— 

la 

h- 

a 

o 

CM 

CM 

H 

fH 






rH 



CO 

VO 

ON 

CO 

CO 

rH 

IA 

CO 

co 

CO 

vo 

c 

0) 

ai 

t- 

ITS 

vc 

LA 

t- 

t- 

CO 

vo 

CO 

vo 

s* 

o 

rH 

On 

ON 

On 

-St 

O 

o" 

CM 

-St 

h~ 


CO 

CM 

00 

H 

CO 

Os 

h- 

-St 

I s - 

rH 

ON 

o 

rH 

CM 

rH 

CM 

rH 





H 




CM 

CM 

rH 

r- 

H 

LA 

rH 

h- 

On 

VO 

CO 


O 

O 

O 

-St 

On 

t- 

CM 

CO 

-st 

vc 

On 

4)1 

• • 

• 

• 

• 

• 

• 

• 

• 

•- 

• 

•• 


CO 

rH 

vc 

h- 

O 

CO 

CM 

VC 

H 

h— 

CO 


IA 

CO 

CO 

vc 

vc 

vc 

On 

CO 

t— 

a 

CO 







rH 

H 

rH 





2 



P* 


00 

On 

CO 

CO 

00 

CM 

CM 

vc 

IA 

ON 

co 

CO 

-st 

rH 

vc 

CM 

ON 

-st 

vc 

H 

LA 


LA 

-St 

CM 

H 

CM 

CM 


h- 

o 


vc 

CO 

vc 

h~ 

vc 

rH . 

H 

hr 

On 

ON 

00 

• 

• 


• 

• 

.*• 

- • 

ON 

LA 

O 

vo 

VO 

-st 

LA 

vc 

t- . 

> 

CM . 

H 

CM 

CM 


U 


4) 













9 

fe 

vo 

CO 

LA 

c 

CM 

h- 

vo 

O 

LA 

N- 

CO 

3 

u 

h- 

o 

O 

H 

CO 

0\ 

-st ' 

LA 

h- 

vc 

vc 

o 

o 

• 

« 

• 

• 

* 

• 

• 

• 

• 

• 

.• 

CO 


o 

CO 

h- 

-st 

On 

H 

CO 

O 

CM 

-st 

CO 


V; 

vc 

rH 

On 

vc 

H 

CO 

CM 

-st 

CM 

co 

CO 

H 

vc 

H 

CM 


rH 

H 

rH 

H 






H 













H 



& 

LA 

•u 

4) 

■P 

On 

CO 

-st 


VC 

45 

h- 

h- 

CO 

e 


U 

• 

• 

• 


CM 

3 

CO 

o 

On 

Vi. 

CM 

-p ' 


CM 

O 

CM 

CM 


rH 

CM 

CM 

CO 


to 





On 

H 

vc 

-st 

h— 

LA 

vc 

-st 

LA 

-st . 

LA 

LA 

• 

• 

■ • 

• 

• 

•/ 

VC 

h- 


CM 

CM 

CO 

H 

O 

vc 

CO 

LA 

-st 

CM 

CM 


















bl 

4) 

4) 


4). 





:L - 

„ •- - • 


« 

•p 

■P 

co 

«P 

-st 

vo 

O. 

CO 

vo 

-it ■ 

C 


a) 


o 

a) 

CM 

; . Ov 

GO 

vo 

vo 

.CO ; ' 

Oi 

O 

c 


• 

E 

• 

• 

• 

• 

• 

‘ ’• 

• 



3 

On 

3 

LA 

o 

o 

CM 

o. 

h- 

H 

•P 

■p 

O 

-p 

rH 

rH 

rH 

CM 

LA 

On 

CO 

-st 

d 

cd 

CM 

«j 

CM 

CM 

CM 

H 



vc 

m 

co 


w 












82 


the sample, and total reflectance is measured in an integrating sphere. 
Percent reflectance is determined as the ratio of energy reflected from 
the sample compared to a standard reflectance material (usually MgO). 

In the field or airborne situation, illumination is more or less hemi- 
spherical, and radiance is measured from a single detector location, 
approximately normal to the panel. At the present time no information is 
available as to the magnitude of differences between laboratory and field 
reflectance measurements, and the assumption was made in this study that 
the differences between the two would not be too great. 

Two different sets of color panels were available for the two 
flight dates. On the July 12 scanner flight date, a set of eight re- 
flectance panels from the University of Michigan's Willow Run Labora- 
tories were deployed at the Purdue University Agronomy Farm. Sample 
patches of each reflectance panel were taken on July 12, and DK-2 spectral 
response in the 0.35 to 2.6 pm wavelength range was determined for each 
sample patch. The set of eight reflectance panels used at the time of 
the July 21 flight was that belonging to LARS. The LARS reflectance 
panels were new at the time of the July 21 flight, and the manufacturer ' s 
curves for DK-2 spectral response in the 0.35 to 2.6 ym wavelength range 
were used. The spectral response curves for the LARS red, green, and 
blue color panels in the 0.35 to 2.0 pm wavelength range as determined by 
the DK-2 spectroreflectometer were plotted (Figure 1 6 ). The DK-2 spectro- 
reflectometer measurements indicate that the gray reflectance panels have 
a relatively constant percent reflectance over the wavelength range from 
0.35 to 2.0 Pm, (Figure 17) whereas the color panels peak sharply in 
certain wavelength regions. 




Figure 16. DK-2 spectra reflectance for red, green, and blue LABS color panels. 





Reflectance Panel 



Wavelength ( /im ) 

spectral reflectance for 5 LARS gray scale panels. 
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Segments corresponding to the eleven wavelength hands of the multi- 
spectral scanner were designated on the DK-2 spectral response curves for 
the reflectance panels. The percent reflectance for each panel was 
determined at the midpoint of the wavelength hand. The percent reflec- 
tance values for the eight reflectance panels for both the Michigan and 
LARS panels were determined and recorded (Table ll). Percent reflectance 
for the color panels was difficult to determine from the DK-2 spectral 
curves in wavelength hands in which spectral response changed notably. 
Therefore, some of the reflectance values for the color panels were 
questionable, and these particular color panels in certain wavelength 
bands were not used in relating scanner response to panel reflectance. 

Evolving Prediction Equations for Relative 
Reflectance from Scanner Data Values 
Plots were made of percent reflectance (Table 11) versus scanner 
data values (Table 10) of the reflectance panels for all eleven wavelength 
bands for both flight dates. Certain panel values were deleted from 
consideration because of saturation of scanner data values or difficulty 
in estimating DK-2 reflectance in certain wavelength bands. Specific 
panels were selected for use in evolving prediction equations for relative 
reflectance from scanner data values (Table 12). 

Percent reflectance (from DK-2 data) versus scanner data values for 
channel 7 of the July 21 flight were plotted (Figure 18). Five panels 
were selected to describe the relationship between percent reflectance 
and scanner data values. Three regression models were theorized: 

1. linear equation through the origin, 2. linear equation with negative Y 
intercept, and 3. second order equation forced through the origin. The 
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Table 11. DK-2 percent reflectance for two sets of color panels in the 
wavelength regions corresponding to channels of the multi- 
spectral scanner. 

Michigan Reflectance Panels - Used on .July 12 Overflight 


Grsy Scale Panels Color Panels 


Channel 

64£ Gray 

32$ Gray 

l6$ Gray 

8$ Gray 

4$ Gray 

Red Green 

Blue 

1 

46 

29 

15 

8 

4 

4 

21 

29 

2 

47 

29 

15 

8 

4 

4 

30 

21 

3 

48 

29 

15 

8 

4 

4 

31 

13 

4 

49 

28 

14 

7 . 

4 

5 

23 

7 

5 

50 

28 

14 

7 

4 

11 

17 

6 

6 

51 

28 

14 

7 

5 

46 

10 

6 

7 

52 

28 

14 

7 

5 

55 

9 

7 

8 

55 

28 

14 

7 

5 

71 

26 

37 

9 

58 

26 

14 

8 

6 

75 

77 

73 

10 

55 

25 

14 

8 

7 

63 

70 

71 

11 

39 

23 

14 

9 

7 

45 

51 

54 


LARS Reflectance 

Panels - Used on July 21 Overflight 





Gray Scale Panels 



Color 

Panels 

Channel 

64$ Gray 

32$ Gray 

16$ Gray 

8$ Gray 

4$ Gray 

Red Green 

Blue 

1 

62 

31 

18 

9 

3 

3 

17 

39 

2 

6 1 

31 

17 

8 

3 

3 

28 

29 

3 

6 1 

30 

17 

8 

3 

3 

34 

16 

4 

6l 

29 

16 

8 

3 

4 

30 

9 

5 

60 

28 

15 

7 

3 

7 

18 

6 

6 

6o 

27 

14 

7 

3 

30 

8 

5 

7 

6o 

27 

14 

7 

3 

62 

Q 

6 

8 

6l 

26 

13 

6 

3 

73 

7 

26 

37 

9 

62 

25 

12 

6 

4 

75 

77 

73 

10 

62 

27 

13 

7 

4 

64 

74 

76 

11 

62 

30 

15 

6 

5 

50 

6o 

65 
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Table 12. Reflectance panels selected for regression analysis on two 
flight dates, 

July 12, 1971 Scanner Flight 


No. of 


Channel 

Observations 



Reflectance 

Panels 

Used 

1 

6 

45? 

85? 

165? 

32% 

Red 

Blue 

2 

6 

h% 

8% 

165? 

32% 

Red 

Blue 

3 

5 

k% 

8% 

165? 

32% 

Red 


4 

6 

k% 

8% 

16% 

32% 

Red 

Green 

5 

6 

k% 

8% 

16% 

32% 

Green 

Blue 

6 

6 

45? 

8% 

165? 

32% 

Green 

Blue 

7 

5 

45? 

8% 

165? 

32% 

Blue 


8 

5 

h% 

8% 

l6% 

32% 

Red 


9 

5 

k% 

8% 

16% 

32% 

Red 


10 

5 

45? 

8% 

165? 

32% 

Red 


11 

5 

k% 

8% 

165? 

32% 

Red 



July 21, 1971 Scanner Flight 


Channel 

No. of 
Observations 



Reflectance 

Panels 

Used 

1 

5 

It 5? 

8% ■ 

165? 

32% 

Red 


2 

U 

45? 

8% 

165? 

Red 



3 

5 

b% 

85? 

165? 

325? 

Red 


4 

5 

' k% 

8% 

165? 

32% 

Red 


5 

6 

k% 

8% 

16 # 

32% 

Green 

Blue 

6 

6 

b% 

85? 

165? 

32% 

Green 

Blue 

7 

5 


85? 

165? 

32% 

Blue 



8 

6 

k% 

8% 

165? 

32% 

645? 

Red 

Green Blue 

9 

8 

k% 

8% 

16% 

32% . 

645? 

Red 

Green Blue 

10 

8 

k% 

85?. 

16% 

32% 

645? 



11 

6 

45? 

8% 

165? 

32% 

645? 




Channel 7 (0.61-0.70 M m) 
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first model had a logical basis in the fact that it included the point 
Where scanner data value equaled percent reflectance which equaled 0, 
but first order regression through the origin failed to provide a satis- 
factory fit to the five data points. The second model provided a 
satisfactory fit to the five data points, but resulted in negative 
estimates of percent reflectance for scanner data values below 31 (of 
which there were two cases among the July 21 Russell plot scanner data 
values). The third model provided a satisfactory fit to the five data 
points, and provided reasonable estimates of percent reflectance for the 
ground cover plots. A second order regression model, forced through the 
origin, was used to describe the relationship between percent reflectance 
and scanner data values for channels 1 through 7 in the visible spectral 
region. 

In the near infrared channels a much different relationship was 
apparent. Percent reflectance versus scanner data values for channel 9 
of the July 21 flight were plotted (Figure 19). A linear regression 
equation best described this relationship, as it did for all the channels 
in the reflective infrared. 

Prediction equations were determined separately for all eleven 
channels for both flight dates (Table 13). These prediction equations 
could be used to determine the relative reflectance, normalized to ground 
reflectance panels , for ground targets in environmental proximity to the 
panels. 
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Table 13. Prediction equations for percent reflectance (relative to 
ground reflectance panels) from scanner data values cali- 
brated to CO and Cl. 

July 12, 1971 Scanner Flight 


y * $i x + e^x 2 


Channel 

fl 

fll 

1 

0.04411 

0.00038 

2 

0.01522 

0.00035 

3 

0.00985 

0.00046 

4 

0.11166 

0.00031 

5 

0.09953 

0.00022 

6 

0.08980 

0.00018 

7 

0.09053 

0.00021 




Channel 

V 

fl. 

8 

- 9.50851 

0.58162 

9 

- 14.67859 

0.50190 

10 

-19.30879 

0.34549 

11 

-9.85590 

0.23758 


July 21, 1971 Scanner Flight 



Y * B-jX + e u x 2 


Channel 

h 

fll 

1 

0.04423 

0.00043 

2 

0.02194 

0.00026 

3 

0.12514 

0.00029 

4 

0.02836 

0.00046 



Table 13* cont 


Channel 

h 

fll 

5 

0.08857 

0.00023 

6 

0.05831 

0.00031 

7 

0.06725 

0.00031 


Y * B + B V X 
O l 


Channel 

fo 

h 

8 

-II.79622 

0.60379 

9 

-16.56261 

1.27868 

10 

-12.58878 

0.76388 

11 

-19.13257 

0.99577 
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Determining Relative Reflectance of Plots Normalized to Panels 
In order to apply the prediction equations (Table 13) to the ground 
cover plots , the scanner data values for the plots had to be obtained 
using C0-C1 calibration. Mew scanner data values vere obtained for the 
plots in the same manner in which they were obtained for the ground 
reflectance panels (Table 14). 

These prediction equations were then applied to each channel of 
data for each of the Russell plots listed in Table 1^. Percent reflec- 
tance values relative to the ground reflectance panels were determined 
for the ground cover plots. Normalized reflectance values were calcu- 
lated and recorded (Table 15). 

Ratio Techniques for Relating Normalized 
Reflectance to Leaf Area Index 

The same ratio techniques described in the previous chapter were 
applied to the normalized reflectance data. The ratios of normalized 
reflectance in channels 8/7 and 9/7 were taken from the values listed in 
Table 15 and related to leaf area index. These ratios and the corres- 
ponding values of LAI were compared (Table 16). Stepwise multiple 
regression analysis was run on the data to determine the relationship of 
ratios of normalized reflectance to LAI for the ground cover plots. 

Results and Discussion 

Evaluation of Normalized Reflectance Curves of Corn Plots 
In order to understand the spectral response from corn canopies it 
is first helpful to get some idea of the individual spectral response of 
corn leaves and the soil background. DK-2 spectral reflectance curves 


Table lU. Mean scanner data values calibrated to CO and Cl for 12 Russell plots for two flight dates. 



Table lU, cont. 
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Table 15. Normalized reflectance (relative to ground reflectance panels) of Russell ground cover 
plots for tvo flight dates. 
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Table 16. Ratios of normalized reflectance in channels 8/7 and 9/7 for 
Russell plots on two flight dates. 

July 12, 1971 Scanner Flight 


ot No. 

LAI 

Ratio 

Channels 8/7 

Ratio 

Channels 

1 

0.01 

1.91 

2.65 

2 

2.89 

11.78 

13.24 

3 

3.48 

15.11 

16.55 

4 

1.30 

6.15 

7.78 

5 

0.01 

2.00 

2.83 

6 

2.34 

10.04 

10.88 

7 

1.11 

6. 77 

8. 06 

8 

0.04 

2.07 

3.44 

9 

3.10 

12.69 

14.08 

10 

1.17 

7.03 

8.84 

11 

2.93 

11.71 

13.08 

12 

3.82 

13.71 

15.09 


July 21, 1971 Scanner Flight 


Plot No. 


LAI 


Ratio 

Channels 8/7 


Ratio 

Channels 9/7 


1 

0.11 

1.65 

2.24 

2 

3.00 

r 


r 

3 

4.28 

Missing 

l 

Missing 


Data 

[ 

Data 

4 

2.31 

t i 


1 i 

5 

0.12 

1.65 

2.17 

6 

3.57 

11.77 

15.28 

7 

2.15 

10.89 

12.73 

8 

0.17 

1.96 

2.61 

9 

3.38 

14.31 

15.09 

10 

2.10 

7.83 

9.37 

11 

3.45 

11.51 

14.66 

12 

4.06 

12.69 

15.94 
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for corn leaves with 80 % moisture content, and two soils in saturated 
and air dry conditions were obtained (Figure 20). The spectral response 
curves for Chalmers silty clay loam, a dark surface soil, and Fincastle 
silt loam, a light surface soil, are shown. Fincastle silt loam is the 
somewhat poorly drained member of the catena of which Russell silt loam 
if the well drained member. The spectral curves for the Russell soil 
should be very similar to those illustrated for the Fincastle soil since 
they have the same surface color and texture and about the same organic 

matter content. It is obvious that moisture content of the soil can 

/ 

greatly affect the spectral response of the soil. The surface soil con- 
dition in a field situation would probably be closer to the spectral 
response illustrated for the air dry soil than that for the saturated 
soil (Hoffer and Johannsen, 1969). 

The curves of scanner data values versus wavelength for three of 
the Russell plots on the July 12 flight date were plotted (Figure 21). 

The plots represent three greatly different ground cover conditions. The 
scanner data values used are the uncalibrated scanner response values from 
the July 12 multispectral scanner mission over the Agronomy Farm (Table 7) 
The wavelength scale is incremented in micrometers on the bottom of the 
graph, with the corresponding midpoints of the channel wavelength bands 
being displayed at the top of the graph. It can be seen that there is no 
relationship between adjacent channels and that the shape of the spectral 
response curves can in no way be related to any familiar response curves 
for green vegetation or bare soil. 
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Figure 21. Unc&libr&ted scanner response curves for 3 Russell plots, July 12. 
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Normalized spectral response curves for three different ground 
cover situations were plotted for two scanner flight dates (Figures 22 
and 23). The normalized response curves of Figure 22 are of the same 
three ground cover plots shown in Figure 21. The same original data were 
used for plotting these curves. The only difference is that the scanner 
data values in the latter have been normalized to relative reflectance, 
using the ground reflectance panels. The relative reflectance values for 
Figures 22 and 23 were obtained from Table 15. The curves in Figure 22 
resemble the DK-2 spectral response curves for green vegetation and bare 
soil. They have the familiar peaks in the green and near infrared wave- 
lengths for green vegetation and the relatively smoothly increasing curve 
for bare soil (Figure 20). The plot with the higher leaf area index has 
a higher response in the 0.72 to 1.3 pm wavelength region and a lower 
response in the O .65 pm chlorophyll absorption region than does the plot 
with a lesser LAI. The response curve in Figure 22 for Russell plot 1, 
with an LAI of 0.01 (essentially bare soil) resembles quite closely the 
response curve in Figure 20 for air dry Fincastle soil. 

In the plots of normalized spectral response curves (Figures 22 and 
23) it is observed that the plots with high percent ground cover have a 
lower response in channels 10 and 11 than plots with lesser ground cover. 
This is probably a result of the spectral response of vegetation from the 
medium ground cover plots being "mixed" with the spectral response of the 
bare soil. 

The normalized response curve for Russell plot 8 (Figure 23) shows 
much higher response throughout the 0.U6 to 2.6 pra wavelength range than 
for Russell plot 1 (Figure 22), even though the ground cover was slightly 



Scanner Channels 



Figure 22. Normalized spectral response curves for 3 Russell plots, July 12. 





Scanner Channels 
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higher on Russell plot 8. Upon further investigation it vas theorized 
that the great differences in spectral response between these two plots 
were not accountable only to ground cover differences. Examination of 
the Hi-Ranger photography taken over these two plots on July 13 and 
July 21 showed that the soil background appeared much lighter on the 
July 21 photography of Russell plot 8. Weather records from the Agronomy 
Farm Weather Station indicated that a long dry period preceded the 
July 21 flight while a rather substantial rain fell the day before the 
July 12 flight. It is likely then, that the great differences observed 
in the spectral response of the low ground cover plots on the two flight 
dates were accountable more to moisture differences than differences in 
ground cover. 

Relative Reflectance Ratios and Leaf Area Index for Two Flight Dates 

The ratios of normalized reflectance in channels 8/7 and 9/7 were 
calculated for the two flight dates (Table 16). These ratios were then 
plotted against leaf area index (Figures 24 and 25). Stepwise multiple 
regression indicated a linear relationship between LAI and both ratios. 
Using the ratio of 9/7 for normalized data, 9 6.4$ of the variation in 

A 

LAI could be explained by the regression equation, Y = -0.7245 + 0.2735X. 
For the ratio of 8/7 for normalized data, 94.1$ of the variation in LAI 
could be explained by the regression equation, Y * -0.5117 + 0.2971X. 

A considerable improvement occurred in the use of the ratio of 
normalized data in channels 8/7 to predict LAI over the use of uncali- 
brated scanner data values in these channels. The procedure of nor- 
malizing the reflectance of the plots to the ground reflectance panels 
apparently was successful in eliminating variations in scanner response 
between flight dates. 
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. /Normalized Reflectance in Channel 9 \ 

0 10 V Normalized Reflectance in Channel 7/ 

Figure 2k. Leaf area index versus the ratio of normalized reflectance 
in channels 9/7 for two flight dates. 
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Figure 25. Leaf area index versus the ratio of normalised reflectance 
in channels 8/7 for two flight dates. 
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Summary and Conclusions 

Multispectral scanner response can be related to the reflectance of 
ground reflectance panels in arriving at prediction equations for relative 
reflectance from scanner data values. This normalization of scanner data 
to ground reflectance panels allows for extension over time of ratio 
techniques for predicting leaf area index. Regression equations can be 
evolved relating leaf area index to the ratios of scanner data values from 
channels 8 and 9 to scanner data values from channel 7. 

Spectral response curves for corn canopies can be determined from 
the derived prediction equations relating panel reflectance to scanner 
response. The spectral response curves for different ground cover plots 
from normalized scanner data show that the various ground cover response 
curves represent a "mixing" of the spectral response from the green 
vegetation and bare soil components. 

The normalized spectral response curves for the ground cover plots 
indicate an increase in reflectance in the 0.72 to 1.3 wm near infrared 
wavelength region with increasing leaf area index. A decrease in reflec- 
tance was observed for the 0.65 um chlorophyll absorption band with in- 
creasing leaf area index. 

Moisture differences apparently had a strong effect on the spectral 
response of the corn canopies on the two flight dates. Soil moisture 
differences greatly affect the spectral response from low ground cover 
plots. 
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The use of ground reflectance panels aids in arriving at normalized 
reflectance values for com canopies. One drawback is the lack of a 
reflectance panel whose reflectance in the visible wavelength region is 
as low as that of a dense com canopy. For this reason, extrapolation 
of data below the known reflectance value of the reflectance panel is 
necessary. This may introduce error in estimating the normalized reflec- 
tance of dense corn canopies in the visible wavelength region. 

The practical implications of using ratio techniques for analysis 
of ground cover sure certain to become apparent in future efforts in 
remote sensing. The orbital perspective of the Earth Resources Techno- 
logy Satellite ( ERTS ) and SKY LAB will provide a general view of agri- 
cultural crops. With the extremely high altitude and coarse resolution 
from space platforms such as these, it is likely that differences in 
vegetative cover will provide the strongest means of discriminating 
between various healthy green agricultural crops. Ratio techniques 
utilizing information from the near infrared and chlorophyll absorp- 
tion regions should prove useful in analyzing relative vegetative 
cover. 
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